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ABSTRACT 

In  this  report  losses  in  systems  of  parallel  round  conductors  are 
studied.  Both  the  normal  skin  effect  loss  and  the  addit  ional  loss  due  to 
the  close  proximity  of  adjacent  conductors  are  considered.  The  results 
obtained  for  the  parallel  conductors  are  used  to  evaluate  the  radiation 
efficiency  of  electrically  small  multiturn  loop  antennas. 


SECTION  I 


ANALYSIS  OF  SYSTEMS  OF  PARALLEL 
ROUND  CONDUCTORS 


I.  Introduction 

In  a  system  of  parallel  conductors  the  distribution  of  current  over 
die  conductor  cress  section  is  determined  by  two  effects — the  normal 
skin  effect  and  a  proximity  effect.  Both  are  the  result  of  die  same 
phenomenon,  eddy  currents  in  the  conductors.  The  former  is  usually 
considered  to  be  the  result  of  the  net  current  in  a  single  conductor  while 
the  latter  is  due  to  die  currents  in  neighboring  conductors.  For  close 
conductor  spacings,  the  distribution  of  current  due  to  the  proximity 
effect  can  cause  an  increase  in  the  ohmic  resistance  which  is  larger 
than  the  skin  effect  resistance  alone,  i.  e.  larger  than  the  ohmic 
resistance  of  the  isolated  conductors. 

The  skin  effect  in  round  conductors  is  discussed  in  most  texts  on 
electromagnetic  theory  [1],  [2],  [3].  The  proximity  effect  has  received 
much  less  attention.  Most  of  the  theoretical  and  experimental  works  on 
the  proximity  effect  deal  with  two  wire  systems  where  the  wires  carry 
equal  currents  in  opposite  directions.  For  examples,  see  the  work  of 
Kennelly  [4],  [5],  Carson  [6],  and  Dwight  [7],  [8],  This  geometry  has  a 
direct  application  in  the  problem  of  wave  propagation  along  parallel  wire 
transmission  lines. 

The  only  investigations  of  the  proximity  effect  in  systems  with 
more  than  two  conductors  appear  to  be  those  done  in  conjunction  with 
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determining  ohmic  resistance  and  Q  of  inductance  coils.  Of  die 
theoretical  treatments,  Butterworth's  discussion  of  die  alternating 
current  resistance  of  cylindrical  conductors  and  solenoidal  coils  is  the  * 

most  thorough  [9],  [10],  [11].  His  work  is  considered  the  standard 

I 

theoretical  approach  and  is  summarized  in  several  places  [12],  [13],  [14]. 

The  experimental  work  of  Medhurst,  however,  indicates  that  Butterworth*  s 
calculations  of  die  radio  frequency  resistance  of  coils  are  not  valid  over 
as  large  a  range  of  parameters  as  expected;  for  certain  dimensions, 
errors  as  large  as  190%  were  observed  [15]. 

In  the  remainder  of  this  chapter,  systems  composed  of  various 
numbers  of  in-line,  parallel  conductors  are  analyzed.  All  the  conductors 
have  the  same  circular  cross  section  and  carry  equal  currents  in  the 
same  direction.  Only  the  high  frequency  case  where  die  currents  are 
confined  to  a  thin  layer  near  the  surface  of  die  wires  is  considered.  This 
report  is  an  extension  of  the  investigation  of  the  two  turn  loop  antenna  * 

reported  in  [16]. 

s 

2.  The  Nature  of  the  Current  Distributions  in  a  System  of  Parallel 
Conductors 

A.  Proximity  and  Skin  Effects 

In  the  system  of  parallel  conductors  illustrated  in  Fig.  1-1  there 
are  two  factors  which  determine  the  distribution  of  current  over  the 
cross  section.  The  first  is  the  normal  skin  effect  which,  for  high 
frequencies,  causes  a  concentration  of  the  current  near  the  outer  surfaces 
of  the  conductor.  This  is  depicted  in  Fig.  l-2a  for  a  single,  isolated. 


FIG.  1-1 


PARALLEL 

SECTION 


WIRES  OF  CIRCULAR  CROSS - 
CARRYING  EQUAL  CURRENTS  IN 


THE  SAME  OIRECTION. 


4 


w 


i 


£ 


f 


% 


i 


> 

j 

t 

i 

!? 

f 

i 

\ 

f 

f 

s 


> 

£- 


<E 

g 

O 

D 

O 

2 

8 


IU 

-J 

e> 


r-  O 

S° 

±  UJ  ^ 
<0  f- 


±  u. 
*  o 

CO 

o  ^ 

te 

g* 

g< 

z  z 


u. 

ll 

UJ 


* 

(0 


CD 

§o 

UJ  8s 

h*  <  Ld 


O  JO 


l 


THE  SHADED  AREA  REPRESENTS  THE  DISTRIBUTION 
OF  CURRENT  ON  THE  SURFACE  OF  THE  WIRE. 

FIG.  1-2 


-5- 


rour.d  conductor.  Secondly,  there  is  an  additional  redistribution  of  the 
current  due  to  the  proximity  effect.  This  is  caused  by  the  magnetic 
field  present  at  any  one  conductor  due  to  the  currents  in  the  other 
conductors  of  the  system.  The  proximity  effect  for  two  parallel,  round 
conductors  carrying  equal  currents  in  the  same  direction  is  illustrated 
in  Fig.  l-2b.  In  the  two  conductors,  the  proximity  effect  forces  the 
current  to  the  outside  edges,  much  as  the  skin  effect  forces  the  current 
to  the  outside  surface  of  the  single  conductor. 


B.  High  Frequency  Approximation  for  the  Current  Distribution  and 

Resistance 

At  sufficiently  high  frequencies  the  skin  depth  d  for  a  good 
conductor  is  a  small  quantity  compared  to  the  cross  sectional  dimensions 
and  most  of  the  current  in  the  conductor  is  confined  to  a  thin  layer  near 
the  surface.  The  magnetic  field  external  to  the  conductor  is  approximately 
the  same  as  the  field  of  a  perfect  conductor  of  the  same  shape  carrying 
an  equivalent  surface  current.  An  expression  for  the  time  average  power 
loss  per  unit  surface  area  of  the  good  conductor,  in  terms  of  the  component 
of  the  magnetic  field  tangent  to  the  surface  of  the  perfect  conductor,  is 


P  =  -»RS  (1^1) 2  Watts/ (meter)2 
c  ^0 

In  terms  of  the  surface  current  Kg  on  the  perfect  conductor 
P  =  -|rS  !  Ks  1 2  Watts/(meter)2 


(1-1) 


(1-2) 


(l-3a) 


where  Rs 


is  the  surface  resistance. 


R 


s 


d  = 


J 


2 


(l-3b) 


If  the  conductor  is  cylindrical  and  Kg  is  an  axial  current  density,  the 
power  loss  per  unit  length  of  the  conductor  is 


P 


!Ks  { 2  dw 


Watts/meter 


(1-4) 


where  the  integral  is  over  the  periphery  of  the  conductor. 

For  the  isolated,  circular,  cylindrical  conductor  of  radius  a 
carrying  total  current  I,  rotational  symmetry  applies.  Equation  (1-4) 
reduces  to  the  familiar  "Rayleigh  formula"  for  the  high  frequency 
resistance  per  unit  length  of  a  circular  conductor 


P  = 


Watts/meter 


(1 -5a) 


tj  ®  i 

“Rayleigh  =  Wi  ‘  2KVTT  Ohms/meter  <l-5b) 

which  is  valid  for 

a/d  »•  1  (1-6) 

9 

With  more  than  one  conductor  present  the  current  distribution  and 
external  fields  for  each  conductor  are  no  longer  rotationally  symmetric; 
therefore,  equations  (1-5)  no  longer  apply.  Further  investigation 
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is  necessary  to  determine  conditions  like  (1-6)  which  insure  that  the 
high  frequency  approximation  expressed  in  (1-4)  is  valid. 

Consider  a  system  of  long,  in-line,  parallel  conductors  carrying 
equal  currents  in  the  same  direction  (Fig.  1-1)  with  parameters  such 
that 

a  «  h  ,  pQa  «  1  (1-7) 

n^c^  «  h^  ,  p^n  c  «  1 

c  >  a  (1-8) 


Neglecting  displacement  currents  as  compared  to  conduction  currents, 

the  axial  component  of  the  volume  current  density  J mz  interior  to  the 
th 

m  cylinder  must  satisfy  the  following  partial  differential  equation  in 
cylindrical  coordinates  (r,  0,  z). 


r 


=  0 


(1-9) 


An  e  1Wt  time  dependence  is  used.  The  solution  to  (1-9),  obtained  by 
the  method  of  separation  of  variables,  which  has  the  desired  symmetry 
and  remains  finite  at  the  origin  is 


J 

mz 


(r,0,z) 


00 

E 

P=0 


M)PC  M  (^fTf-) 

mp  p  *  d 


i0p^| 


') 

8  cos(P0) 
(1-10) 


where  M  and  0^  are  the  modulus  and  phase  of  the  Kelvin  functions 

p  p  r 

(ber  +  i  bei  )  [17,  p.  379].  The  C  are  functions  of  z  only. 

P  P  P 


The  total 
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current  at  a  given  cross  section  of  the  conductor  is  therefore 


m 


r  r  2fads  cz  *’  ei^f) 

w  *  f  /  =  ~rT  CmOMlW?  3->  « 

r=0  0=-jt  *2  8 


(1-11) 


and 


I  (*) 

m 


'mO 


e* 


^/?irads 


(1-12) 


The  volume  density  of  current  extrapolated  to  the  surface  of  the 
conductor  is 


J  (a  d  z)  = 
mz  *  ’ 


_  8  e  s  S 


00 


l1  +  ^  anip  cos  (p0)] 


(1-13) 


P  1 


With  (1-12)  and  (1-13)  substituted  in  (1-10),  the  current  density  becomes 

'm<*>  M0 I^T21  -‘I  ®o('^'tA)  - 


J  (r,  0,  z)  =  — - — 

mz  ^7radgM1(>/?|-) 


oo 


1  +  )  a’  - - - - S —  cos  (pQ) 

I  U  mp 

I  S  r  8 

•>l  9  0(V?F>  -  ®0^f£')| 


(1-14) 
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Where  •  =  a  -  r  is  the  radial  distance  into  the  conductor  from  the  snrfice. 
In  the  present  analysis,  the  coefficients  a^  are  assumed  to  be  complex 
numbers. 

When  the  current  distribution  at  file  surface  of  the  conductor  is 
sufficiently  smooth,  a  finite  number  q  of  the  Fourier  series  terms  in 
equation  (1-14)  are  adequate  to  approximate  the  current  density.  If,  in 
addition,  file  frequency  and  conductivity  are  high  (a/d#  »  1,  r  «  a)  file 
large  argument  asymptotic  formulas  for  and  8^  apply  [l?].  In¬ 
serting  these  into  equation  ( i -14)  yields 


(1-15) 


which  simplifies  to 


J 

mz 


(r,0z)  = 


(1-16) 


(1-17) 
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Witfc  these  conditions  satisfied,  file  current,  although  non-uniform  in  B , 

is  confided  to  a  thin  layer  near  the  surface  much  as  in  file  case  of  an 

tJl 

isolated  cylinder.  The  power  dissipated  per  unit  length  in  the  ir» 
conductor  is  then 


P* 

m 


7  f  f 

r=0  d-t 


rdrdfl 


l*(z)  R* 

m 


4*a 


M  E' 

pi 


I2l 

mp‘  1 


Watts/meter  (1-18) 


If  the  cylinders  are  now  made  perfectly  conducting,  the  current 
th 


on  the  m  cylinder  will  be  of  the  form. 

I  (*> 


I  (z) 

K  (0,z)  --  ~ —  * 
mz  ’  2ya  m 


«„(») •  -gr  'i +  Yu  *mp co* (p9)| 


q 

p=i 


(1-19) 


where  gm(0)  is  the  normalized  surface  current  density.  Using  the 


approximation  expressed  in  equation  (1-4),  the  power  loss  per  unit  length 

for  a  good  conductor  expressed  in  terms  of  the  coefficients  a  for  the 

mp 

perfect  conductor  is 


i  I  (z), 

P  -*  I  /  m  \ 2 
m  "  2  (  2jra  }  R 


g^(0)ad  9 


iL{z)r8 

TCI 

4jra 
Watts/meter 


[•♦iXX/1 

P=i 


(1-20) 
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For  large  values  of  a/df  ftii  expression  is  a  good  approximation  to  the 
correct  relation^  equation  (1-18),  that  is 


a 

mp 


amp 

P» 

m 


provided 


Pd,  2 

<-jr> 

- 7J- 

a  -  (f»V> 


«  1 


a/d#  »  1 


P  =  1,  2--q 


(1-21) 


(1-22) 


The  first  term  in  equation  (1-20)  is  the  power  loss  in  the  m4*1 

conductor  due  to  the  net  current  I  in  that  wire.  This  is  the  normal 

m 

skin  effect  loss.  The  sum  in  (1-20)  represents  the  loss  due  to  nonuniform 

currents  induced  by  other  wires  in  the  system.  J.t  is  the  additional  loss  in 

the  m4^  wire  due  to  the  proximity  effect.  Since  the  coefficients  a^  in  the 

sum  are  a  function  of  the  net  currents  Ij  in  all  wires  of  the  system,  the 

1  2 

equation  for  P__  cannot  be  written  as  P  =  R  (-y  I  )  if  R  is  to  be  only 
m  m  m  m  m 

a  function  of  the  physical  parameters  of  the  system.  As  a  result,  the 

usual  circuit  definition  of  the  ohmic  resistance  of  each  wire  (R  - 

'  m 

1  2 

P  /-yl  )  makes  no  sense, 
m  c  m 

When  all  conductors  carry  the.  same  total  current  at  each  cross 
section  the  ohmic  resistance  per  unit  length  of  the  system  of  wires  is  a 
useful  quantity.  Using  the  series  definition  of  the  current  (1-19)  the 
ohmic  resistance  per  unit  length  of  a  system  or  n  parallel  wires  is  given 
by 
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8  “  q 

R  =  52  11  4  £  ]C  ^mp^3  ohm*/meter  (1-23) 

m=l  p=l 

If  the  separation  between  conductors  is  large  enough  that  each  can  be 
considered  as  isolated  from  the  others,  (1-23)  becomes 

Ro  s  E"  =  "  ^Rayleigh  Ohn«/m.ter  (1-24) 

The  additional  ohmic  resistance  per  unit  length  due  to  the  proximity 
effect  is  then 

n  q 

Rp  =  R  -  Rq  =  ^  lampi^  Ohms/meter  (1-25) 

m=  1  p- 1 

Normalized  quantities  are  useful  when  comparing  different  configurations 
of  conductors. 


W  ■  Z|I+2  Z^mpl2)  (1-26' 

m  - 1  p-- 1 


R 

_P 

R 


1 

2 


t  i 

m=  1  p=  1 


(1-27) 


In  the  present  analysis  a  smooth  conductor  with  a  uniform  surface 


resistance  is  assumed.  Recent  research  by  A.  Sanderson  [18]  indicates 
that  surface  roughness  in  the  form  of  scratches  transverse  to  the 
direction  of  the  current  can  significantly  alter  the  equivalent  surface 


resistance.  The  type  of  wire  used  in  practical  applications  is  usually 
formed  by  a  drawing  process,  such  as  drawn  copper  wire.  Surface 
scratches  are  in  the  same  direction  as  the  current  flow  and  are  expected 
to  increase  the  ohmic  loss  much  less  than  equivalent  transverse 
imperfections  would.  Calculations  using  Sanderson's  theory  indicates 
that  surface  roughness  can  be  ignored  at  the  frequencies  of  interest 
(<  100  MHz. )  in  this  study. 

3.  Formulation  of  the  Integral  Equations  for  the  Transverse  Current 

Distributions 

Consider  each  of  the  long,  parallel  cylinders  in  Fig  1-1  as  being 
perfectly  conducting.  The  surface  current  density  on  the  l**1  conductor 
is  then 

Kl(0,’z,)  =  2^^')  *<*’)  ,  f  =  1,  z  -  -  n  (1-28) 

The  dimensionless  quantity  gj(0')  is  the  normalized  surface  current 
density.  In  (1-28)  the  same  z'  dependence  f(z')  is  assumed  for  the 
current  distributions  on  all  cylinders.  The  conductors  are  composed  of 
three  sections;  the  length  z-d^z'^z+d  and  the  two  end  sections 
z+d^z'<h,  -h<z'  <z-d.  In  addition  to  the  inequalities  presented  in  (1-7) 
and  (1-8)  the  following  constraints  are  placed  on  the  length  d 

PQd«l  (l-30a) 

d2»n2c2  (l -30b) 

This  makes  the  current  distribudons  at  every  cross  section  along  the 
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length  z-d^z*  <z+d  approximately  die  same. 

Kf(0»,z‘)  =  ^Igl(a,>  f(z)  »  *-**»' Orfd  0-31) 


The  Helmholtz  integral  for  the  vector  potential  component 

*♦  ^ 

A  (r ,0,  z)  at  a  point  just  off  the  surface  of  the  m  conductor  is 
z 


0-32) 


where 

Rmf  r  {(z-z')Z  +  r2^]1/2  [(z-z')2  +  4(m-l)2  c2  + 

r2  +  a2  -  2ar  cos  (0  -0')  +  4(m-f)c(r  cos0  -  acosS')]1^2 

0-33) 

If  terms  of  order  or  less  are  neglected  in  the  first  integral  and 
setting 

Rmf  1  [<z-z')2  +  4(m-f)2  c‘’]l/2  (1-34) 

in  the  last  two  integrals,  equation  (1-33)  reduces  to 
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v» » 5  «  J  /  e[Ci d9,d" 

z'?z-d  a»=-w  i=1 

/  z?*d  ^  \  n  r  *ok«’>* + 4<m"1^2 

I  /  /  i2^f(*')  \  '  I  i  2  i  1j 

\  J  J  I  Tl  Kz*z'l  +  4(m-l)Z  c  ]ly 

\z'=-h  z'= z+d/  L 


(1-35) 


The  z*  integration  in  the  first  integral  can  be  evaluated  directly  [19, 
p.  50,  200.01].  The  result  is 


/  ^ 

z'=z-d 


—  dz*  =  2  sin  h"1  (-r-)  =  2[!n(r  .)  -  in(2d)]  (1-36) 

ml  ml 


2  2  2 

where  terms  of  order  n  c  /d  are  dropped  in  the  last  expression. 
With  (1-36),  equation  (1-35)  becomes 


v".9.'' 5  A  -2«*> 

8;r 


+  4jt  nln(2d)  +  A^fz) 


/&*■ 

. .  !-  1 

0’  =  -t>  1 


)ln(rm#)]de' 


(1-37) 


The  term  AJ^tz)  represents  the  last  two  integrals  in  (1-35). 

The  normalized  surface  current  density  is  given  by  the  boundary 
condition 


-Zita.  8  A  (r,0,z) 

gm(0)  =  -  —£5 - 

HQIf(z)  9  r 


(1-38) 


where 


Comparing  this  with  the  following  definition  of  the  Dirac  delta  function 


5(0)  =  -  lim  [—A  -z  ] 

v  A~*Q  A2+0Z 


(1-44) 


it  is  evident  that  in  the  limit  the  integrand  becomes 


Integrand  =  [ff  5(0)  + 


(1-45) 


Substituting  (1-45)  into  (1-39)  and  rearranging  yields 


«m|S) 


•f  n 

9’=-ir  fcV 


(1-46) 


where 


K  ,(0,0») 
m,  i ’ 


1  t  2(m-l)  c/a  coa  0  -  cos  (0-01] 

—5 


(1-47) 


Symmetry  about  the  center  of  the  system  of  wires  requires  gn+j_m(0) 
gm(tr-0)  which  reduces  the  number  of  terms  in  (1-47)  to  n/2  for  n  even 
or  (n+l)/2  for  n  odd. 


n  even 


gm(0) 

m 


K  n4,  w(0,7r-0»)  g  (0')d0'  + 
m.n+i-m  __ 

'  IXl 


n/2 


0^-ir 


ir. 


0  -  -7T 


^Km,i(0’0,)  +  Vn•M-f(e'7^"e,)lgf(0,)d9,  ,  -  =  1,  2-n/2 


(1-48) 
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n  odd 


g(0)  = 
m 


y  y 

I  K  ^  (0,9-0*)  g  (d*)dd»  1+-  J  / 

I  m,  n+l-m  ’  7  6m  7  ff  f  '  / 

»*  i-i  ^ 

»'=-»  Vm  9’=-" 


(1-49) 


where 


h{m)  =1  m  |  (n+l)/2 

=  0  m  =  (n+l)/2 


(1-50) 


Equations  (1-48)  and  (1-49)  represent,  respectively,  a  system  of  n/2  and 
(n+l)/2  coupled  integral  equations  whose  solutions  are  the  desired 
surface  current  densities,  gm($)  • 


4.  Solution  for  Two  Conductors 

The  simplest  geometry  for  examining  the  proximity  effect  is  two 
parallel  circular  conductors  carrying  equal  currents.  Exact  expressions 
for  the  current  distribution  and  ohmic  resistance  for  this  simple  case 
are  given  in  Technical  Report  No.  612  [16],  The  normalized  current 
distribution  on  the  two  wire  system  is  graphed  in  Fig.  1-3  for  variois 
conductor  spacings  c/a.  The  distributions  will  prove  useful  in  developing 
approximate  solutions  for  systems  with  two  or  more  conductors. 
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5.  Approximate  Solution  of  the  Integral  Equations  for  Two  or  More 
Conductors 

A.  The  Method  of  Undetermined  Coefficients 

1.  Reduction  to  a  Set  of  Algebraic  Equations 

For  systems  with  more  than  two  circular  conductors  approximate 
methods  must  be  used  to  obtain  the  current  distribution  and  resulting 
ohmic  resistance  per  unit  length  of  the  system.  In  this  section  one  such 
method,  undetermined  coefficients,  is  applied  to  the  previously  derived 
system  of  integral  equations  (1-48),  (1-49). 

As  the  previous  analysis  suggests,  a  trigonometric  series  is  the 
natural  choice  for  an  expansion  to  represent  the  normalized  surface 
current  density. 


q 

gm(0)  1  +  Yh  ampCOS  (P<))  (,_5l) 

p^l 


Further  evidence  for  this  selection  is  found  by  examining  the  exact 

solution  for  the  two  wire  case.  A  Fourier  analysis  of  the  current 

distribution  (Fig.  1-3),  for  the  limiting  case  c/a  -  1,  indicates  that  the 

first  two  cosine  terms  in  the  series  are  adequate  to  predict  the  correct 

value  of  the  additional  loss  due  to  proximity  R  /R  to  within  1%.  For 

p  o 

large  spacings,  c/a  »  1,  the  current  distribution  is  of  the  form  1  +  a^  cos  0 
as  is  evident  from  Fig.  1-3.  This  last  statement  is  also  true  for 

systems  with  more  than  two  wires  and  is  easily  understood  if  the 
magnetic  field  due  to  external  currents  is  considered  a  constant  over  the 
cross  section  of  each  conductor.  The  magnetic  field  normal  to  the 
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axis  of  tiie  m^1  conductor  in  a  system  of  n  conductors  would  be 


Hq1 

J  =  -- . 
my  4ffc 


** 

H  (n^T) 


(1-52) 


i=l 

f/m 


and  the  resulting  current  distribution  becomes 


5  1  +  <?  E  (Sb) 1  COS  9 
1=1 
f/m 


(1-53) 


Substituting  (1-51)  into  the  integral  equations  (1-48),  (1-49),  one 
obtains 


n  even 


2  “mp  I"  008  'P®'  + 

P=1 


I  K 

7T  /  m, 
0'  =  -ff 


n+l-m*0'0'*  cos  (p0,)  d0'J 


n/2  q  |  V 

L  IJ%|f  J  !Km,i<«.«'>  +  <-1>PVn+i.1(e,9')]co.(pe')de 

l/m  P'  e'"-n 


i/j 


Fm,n+l-m<e-e'l  +  2^  [Km,i<e-P’)  +  Vn+>-i(e’e’>lide’ 

1=1 
f/m 


m 


=  1,  2,  —  n/2 


(1-54) 
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n  odd 


y  a  [ 
Z,  mpl 

p^l 


-cos 


p0  m/(n+l)/2  ,  nP  / 

.y  JLiL  h(m)  /  K  . 
2p0  m=(n+l)/2  7  J  ’  1 


0')cos(p0')dtf '] 


0'=-ir 


(n-l)/2  q  j  r 

1  Sa<P  »  /  !KtnIJ(e’9')+("1)i>Km>n+l-l(9’e')'coslpe,,de' 

f/m 

q 

2  Vl  -  ^  (n+1)/?.^*  ^')cOS^^,)d^'^ 

P  “  P  0<U 


I  /  |h(m)K  (fl, 

nil  m,  n+l-m  ’ 

0'=-7T 


(n-l)/2 


J 


*  2 

t  l  ' 


t/ 


d0  m  1,  2, - (n+l)/2 

(1-55) 


m 


where  the  same  number  of  harmonic  terms  q  is  used  to  represent  the 
surface  current  on  all  conductors.  Due  to  symmetry  about  9  -  ff/2,  only 
even  harmonics  appear  on  the  center  conductor  of  a  system  with  an  odd 
number  of  conductors. 

The  definite  integrals  in  equation  (1-54)  and  (1-55)  are  of  the 


form 


I((l,m-f,p)  -  - 


[l+2(m-l)(c/a)cos  9  -  (cos  9  cos  9' 
[4(m-?)2(c/a)2  +  2  +  4(m-f)(c/a)cos  0 


t  sin  U  sin  Q')\  cos(pd')  dO1  _ 

-  (4(m-!)(c/a)  +  £'cos(0»))  co.-  O'  -  2  sin  9  sin  f V 


l 1  -5ol 
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Appendix  A  contains  a  detailed  evaluation  of  this  integral,  the 
which  are 


1(0,  m-f,  p) 


~  2~  ,p»l  IA»2  +  Bs  +  C]  ,  p  =  1,  2,  -q 

(1-s  )(-s)y 

- -y  [Bs  +  C]  ,  p  =  0 

s(l-s*) 


where 

s  =  (4(m -1)^  (c/a)^  +  1  +  4(m-i)  (c/a)  cos  0)*^2 
A  =  cos  (0  -  (p  -  1)  i// ) 


B  =  2(1  +  2(m-£)  (c/a)  cos  0)  cos  (pi//) 
C  =  cos  (0  +  (p  +  1)  \f/  ) 


"  tan  (2(m-l)(c/a)  +  cos  0  = 

xf/  =J 

(tan_1 

The  principle  value  of  tan  *  is  used  in  (l-57f).  For  the  case  n 

system  of  equations  (1-54)  with  (1-57)  becomes 

<J  n/2  q 

2  amp  I"  cos  (p0)  +  (_1)P  1  2m-n-i,p)]  +  ^  afp 

P= 1  f  =  1  p=  1 

f/m 


[1(0,  m-f,  p)  +  (-1)^  1(0,  m+f-n-1,  p)j  =  -  {  1(0,  2m-n~l,  0) 


+ 


n/2 

I 

jUI 

f/m 


[1(0,  m-f,  0)  +  1(0,  m+f-n-1,  0)] } 


results  of 

(l-57a) 

(l-57b) 
(1  -57c) 
(l-57d) 
(l-57e) 

(1-57  f) 

even  the 


m  -  1,  2,  --n/2 


(1-58) 
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I 


I 


I 


[ 

; 


* 


This  is  a  set  of  n/2  equations,  one  for  each  conductor,  involving  qn/2 

unknowns  (a  )  the  coefficients  of  which  are  functions  of  the  variable  6. 
mp 

In  order  to  solve  the  system  a  set  of  qn/2  conditions  is  necessary.  Two 
procedures  which  yield  such  conditions  were  used,  the  methods  of 
collocation  and  least  squares  [25]. 


2.  Solution  by  the  Methods  of  Collocation  and  Least  Squares 
In  the  method  of  collocation,  the  am^  are  chosen  so  that  equation 
(1-77)  is  satisfied  exactly  at  q  points  (0  ^  TT,  k  1,2,  --q)  on 

each  conductor,  more  specifically  for  n  even 

q 

2  amp  I"  ““'P'W  +  <-1>P  I(9mk'  2m-n-l,p>l 

p.l 

n  q 

1=  1  P-  1 
i/m 

n/2 

+  2  P<emk-n’-,-0>  +  IWmk-m+|-”-i-0») 
i^l 
i/m 


m  =  1,2,  --n/2 

k  -  1,2,  --q  (1-59) 


With  the  definition  of  the  new  variables 
( 1  -59)  becomes 


.ml 

V 


and  s 


mk* 


equation 


q 

V  .m 

Z  amp  kp 

P'l 


n/2 

2 

f-i 

i/m 


q 


2 


mk 


il-oO) 
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In  matrix  form,  the  system  of  algebraic  equations  (1-59)  is  now 


T11  T12  Tln/2 

- 1 

> 

M 

_ 1 

S1  “ 

T21  T22  T2n/2 

A2 

S2 

Tn  .  Tn  7  '  "  Tn/2  n/2" 

V 

3 

\ 

CN> 

1 _ 

-Sn/2" 

(1-61) 


where 


'‘ii 

£ 

l\z 

•  '  -  *lq 

- 1 

tij 

12 

1 

1 

1 

T..  = 

li 

*21 

*22 

-  -  ‘  4q 

,  T..  = 

tij 

*21 

tM 

22 

cr 

.^rvj  . 

i  ■ 

1  . 

i 

t1  ? 
q2 

1 

1 

1 

IT*- 

•Q 

1 _ 

tij 

Lql 

tij 

q2 

-  -  -  J 

qqJ 

(1  -62a) 


(l-62b) 


r  i 

“  • 

an 

8il 

ai2 

si2 

« 

si- 

• 

aiq 

s. 

iq 

■  m 

m  m 

(l-62c) 

(l-62d) 


For  the  case  n  odd,  a  similar  matrix  equation  results.  All  the  elements 
of  the  S  and  T  matrices  are  real;  therefore,  equation  (1-61)  can  be 


written  as  two  separate  equations 
IT]  [Arj  =  [S] 


(l-63a) 
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[T ]  [A1]  =  0  (l-63b) 

t  i 

where  A  contains  only  the  real  part  of  the  coefficients  Re(amp)  and  A 

contains  only  the  imaginary  part,  lm(amp).  For  a  unique  solution  of 

(l-63a),  the  T  matrix  must  be  nonsingular.  A  nonsingular  T  matrix 

indicates  a  trivial  solution  for  the  homogeneous  equation  ( 1  -63b),  i.  e. 

all  Im(a  )  -  0.  The  normalized  current  distributions  g  (0)  are 
mp  6m 

therefore  real  quantities. 

The  method  of  least  squares  differs  from  collocation  in  that  the 
amp  are  chosen  in  such  a  manner  that  equation  (1-58)  is  satisfied  in  a 
least  squares  sense  over  the  interval  0  <0  <  jr  rather  than  satisfied 
exactly  at  specific  points,  namely 


/*(  q  q 
J  I  2  amp  l‘CO8<P0)  +  (-1)P  1(0.  2m-n-l,p)]  +  £  ^  alp  IW,  m-f,  p) 

oJ  p-1  i;1  p-1 

f*m 


n/2 


+  (-l)p  1(0,  m+f-n-1,  p)]  +  1(0,  2m-n-l,  0)  +  ^  [1(0,  m-f,  0)  +  1(0,  m+f-n-1,  0)  [ 


f=l 

f/m 


d0 


-  minimum  m 


=  1,2  ---  n/2 


(1-64) 


Differentiating  the  left  hand  side  with  respect  to  each  coefficient  a 
and  setting  the  results  equal  to  zero  yields 


mp 


J  [-cos(k0)  +  (-l)kI(0,2m-n-l,k)j 
0-0 

n/2  q 


Y  a  [-cos(p0)  +  (-1  )PI(0,  2m-n- 1,  p)  1 
/ .  mp1 

P  1 


t  ^  aip[I(0,m-f,p)  +  (-1)1’  1(0,  m+f-n-1, p)j  +  1(0,  2m-n-l,  0) 


II  p=  1 
f/m 
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n/2 


1=  1 
f/m 


d0  =  0 


+  ^  [l(0,m-i,  0)  +  1(0,  m+i-n-1,  0)] 

m  =  1,  2,  —n/q 
k  =  1,2,— q 


(1-65) 


After  rearranging  terms  and  performing  integrations,  (1-65)  becomes 

q  i  7r 

j  amp  J  2  ®  (k>  P)  "  /  [( ~  1  )P  cos (k0) 1(0,  2m-n-l,  p)  +  (-l)k  cos(p0)  1(0,  2m-n-l,  k) 
p=l  I  0=0 

n/2  q 

-  (-l)p+kI(0,  2m-n-l,p)  1(0,  2m-n-l,  k)]  d 6 


1  1 


Z,  aip 


f=l  P=1 

f/m 


J  [cos(k0)  -  (-l)kI(0,  2m-n-l,  k)][I(0,  m-f,  p)  +  (-1)PI(0,  m+f-n-1,  p)]d0 
0=0 


J  [cos(k0)  -  (-l)kI(0,  2m-n-l,  k)] 
0=0 


n/2 


1(0,  2m-n-i,  0)  +  V  [1(0,  m-f,  0) 


f=l 

f/m 


+  1(0,  m+f-n-1,  0)] 


d0 


(1-66) 


m  =  1,2, - n/2 

k  =  1,2,  —  q 


which  can  be  written  as 


q 

I 

p=l 


,m  . 
a  t.  + 
mp  kp 


n/2  q 

I  I  a-  ‘ 

f=l  P=1 

f/m 


mi 
fp  "kp 


=  s 


mk 


(1-67) 
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The  variables  t^,  and  enter  the  matrix  equation  (1-62)  in 

the  same  manner  as  in  the  method  of  collocation. 


B.  Numerical  Results 

1.  Comparison  of  the  Two  Methods  of  Solution 
For  the  collocation  solution,  the  same  matching  points 
=  0^,  0  ^  0^  <  n)  were  used  on  all  cylinders  except  the  center 
cylinder  in  a  system  with  n  odd.  The  current  on  the  center  cylinder  has 
symmetry  about  0  =  ir/2,  n;  therefore  only  points  in  the  first  quadrant 
are  needed.  These  were  chosen  to  be  0^/2.  Several  different 
combinations  of  matching  points  were  used  in  (1-59)  and  the  resulting 


matrix  equation  (1-61)  was  solved  for  the  coefficients  a^  using  a 
standard  Gaussian  elimination  algorithm  [26].  The  additional  ohmic 
resistance  due  to  the  proximity  effect  Rp/RQ  was  calculated  from  (1-27) 
for  various  numbers  of  harmonic  terms  q.  No  particular  distribution 


of  points  gave  a  best  rate  of  convergence  of  Rp/RQ  f°r  *11  numbers  of 
conductors  and  spacings.  The  final  set  of  matching  points  settled  on  is 


For  an  even  number  of  harmonics,  the  points  are  equally  spaced  and 
internal  to  the  region  0  <  0^  <  n.  With  an  odd  number  of  harmonics,  a 
slightly  better  rate  of  convergence  was  found  when  the  set  of  matching 
points  did  not  include  0^  --  tt/2. 
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The  least  squares  procedure  requires  the  evaluation  ot  the  de¬ 
finite  integrals  in  equation  (1-66).  Due  to  the  complexity  of  the 
1(0,  m-f,  p)  functions  in  the  integrand,  an  exact  evaluation  wai-  unob¬ 
tainable  and  approximate  numerical  integration  necessary.  A  typical 
integration  from  (1-66)  was  performed  using  three  different  numerical 
integration  routines:  Romberg,  Simpson's  rule,  and  Gauss  quadrature. 
The  six-point  Gauss  quadrature  formula  [27]  required  the  least  time  for 
the  desired  accuracy.  The  interval  0  <0  ^  y  was  divided  into  k+1  or 
p+1  panels,  whichever  was  larger,  and  the  six-point  formula  applied  to 
each  panel.  With  the  integrals  evaluated,  the  resulting  matrix  equation 
(1-61)  was  solved  using  the  same  algorithm  as  for  the  collocation 
solution. 

A  comparison  of  the  two  methods  is  presented  in  Fig.  1-4,  where 

Rp/RQ  and  the  computation  time  for  the  1.  B.  M.  360/65  computer  are 

graphed  as  a  function  of  the  number  of  harmonic  terms  used  in  the 

solution.  The  results  are  for  4  cylinders  with  a  spacing  c/a  =  1. 10. 

Least  squares  is  the  more  elegant  of  the  two  procedures,  converging  to 

the  limiting  value  of  R^/Kq  when  the  number  of  harmonic  terms  is  less 

than  half  that  required  in  the  collocation  solution.  In  terms  of 

computation  time  the  collocation  method  is  much  faster — roughly  6q 

times  faster  for  a  given  number  of  harmonic  terms.  Thus  the  limiting 

value  of  R  /R  is  obtained  in  about  one -tenth  the  computation  time 
p  o 

needed  for  the  least  squares  solution.  Similar  time  savings  are  found 
for  other  numbers  of  cylinders  and  spacings.  For  this  reason,  the 
majority  of  the  calculations  for  this  work  were  do*ie  by  the  method  of 


collocation.  Listings  of  the  computer  programs  for  both  methods  are 
in  Appendix  B. 


2.  Transverse  Current  Distributions 

The  number  of  harmonic  terms  used  for  the  current  distribution 

on  a  given  system  of  conductors  was  determined  by  observing  R^/Rq. 

If  increasing  the  number  by  two  produced  less  than  a  0.  10%  change  in 

R  /R  ,  the  number  of  terms  was  deemed  sufficient.  The  normalized 
p  o 

surface  current  densities  g(d)  for  systems  with  3,  4,  5  and  6  conductors 
and  various  spacings  c/a  are  plotted  in  Figs.  1-5  through  1-8.  The 
distributions  for  2  conductors  are  not  plotted,  since  they  are  identical 
to  those  in  Fig.  1-3.  In  systems  with  three  or  more  closely  spaced 
cylinders  there  are  both  positive  ar.d  negative  currents  on  the  surface 
of  the  outer  conductors.  These  currents  in  opposite  directions  add 
nothing  to  the  net  current  in  the  wire;  they  just  increase  the  ohmic  loss. 
When  the  spacing  between  cylinders  is  very  close,  currents  on  adjacent 
surfaces  of  two  conductors  tend  toward  equal  values  with  opposite  sign; 
for  example:  for  4  wires,  spacing  c/a  =  1.1,  g (it)  =  -  2  on  cylinder  1, 
while  on  cylinder  2,  g(ir)  =  +2. 


3.  The  Additional  Ohmic  Resistance  Per  Unit  Length  Due  to  the 
Proximity  Effect 

Computed  values  of  the  additional  ohmic  resistance  per  unit  length 
due  to  the  proximity  effect  R^/Rq  for  systems  with  various  spacings  c/a 
and  up  to  eight  conductors  are  presented  in  Fig.  1-9  and  Table  1-1. 
Calculations  of  R^/Rq  were  not  made  for  extremely  close  spacings,  i.  e. 


FIG.  1-6  NORMALIZED  SURFACE  CURRENT  DISTRIBUTION  FOR 
FOUR  WIRES 
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Table  1-1.  Normalized  Additional  Ohmic  Resistance  Per  Unit  Length 
Due  to  the  Proximity  Effect  R^/Rq. 


Spacing  c/a 

Number  of  Conductors 

2 

3 

4 

5 

6 

7 

8 

1.00 

0.  333 

1.05 

0.316 

0.  748 

1.231 

1. 10 

0.  299 

0.643 

0.  996 

1.347 

1.689 

2.  020 

2.  340 

1.  15 

0.  284 

0.  580 

0.  868 

1.  142 

1.400 

1.693 

1.872 

1.  20 

0.  268 

0.  531 

0.  777 

1.002 

1.210 

1.401 

1.577 

1.25 

0.  254 

0. 491 

0.  704 

0.  896 

1.068 

1.224 

1.365 

1.30 

0.  240 

0.455 

0.644 

0.  809 

0.  956 

1.086 

1.203 

1.40 

0.214 

0.  395 

0.546 

0.  674 

0.  784 

0.880 

0.  965 

1.50 

0.  191 

0.346 

0.470 

0.  572 

0.  658 

0.  732 

0.  796 

1.60 

0.  173 

0.  305 

0.  408 

0.  492 

0.  561 

0.620 

0.  670 

1.  70 

0.  155 

0.270 

0.353 

0. 428 

0. 485 

0.532 

0.573 

1.  80 

0.  141 

0.241 

0.316 

0.375 

0.  423 

0.462 

0.  495 

1.90 

0.  128 

0.216 

0.  281 

0.  332 

0.372 

0. 405 

0. 433 

2.  00 

0.  116 

0. 195 

0.  252 

0.  295 

0.  330 

0.  358 

0.  392 

2,  20 

0.  098 

0. 161 

0.  205 

0.  239 

0.  265 

0.2S6 

0.  304 

2.  40 

0.032 

0.  135 

0.  170 

0.  197 

0.217 

0.  234 

0.247 

2.  50 

0.077 

0.  124 

0.  156 

0. 180 

0.  198 

0.213 

0.  225 

2.  60 

0.071 

0.  114 

0.  144 

0. 165 

0.  182 

0.  195 

0.  206 

2.  80 

0.  061 

0.  098 

0.  123 

0.  141 

0.  154 

0. 165 

0.  174 

3.  00 

0.  054 

0.  085 

0.  106 

0.  121 

0.  133 

0.  142 

0.  150 

3.  50 

0.  040 

0.  062 

0.077 

0.  087 

0.095 

0.  101 

0.  106 

4.  00 

0.031 

0.043 

0.  058 

0.066 

0.072 

0.076 

0.  080 

FIG.  1-9  THE  ADDITIONAL  OHMIC  RESISTANCE  PER  UNIT  LENGTH  OF  A  SYSTEM 
OF  PARALLEL  WIRES  DUE  TO  THE  PROXIMITY  EFFECT. 
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3  and  4  wires,  c/a  less  than  1.  05;  5  or  more  wires,  c/a  less  than  1. 10> 
The  reason  for  this  will  be  evident  after  a  closer  examination  of  the 
approximation  already  made. 

In  the  limit  as  c/a  approaches  1.  0,  the  surface  current  develops 
large  spikes  at  adjacent  points  on  successive  cylinders.  This  is 
illustrated  for  3  wires  with  s pacings  c/a  =  1. 10,  1. 05,  and  1. 01  in  Fig. 
1-10.  For  wires  with  finite  conductivity,  a  change  in  the  form  of  the 
current  distribution  in  file  radial  direction  is  expected  to  accompany 
these  areas  of  high  current  density.  As  a  result,  the  radial  decay  rate 
will  differ  from  the  high  frequency  skin  depth  d  in  these  regions.  This 
is  basically  the  same  idea  expressed  in  equation  (1-22).  Spikes  in  the 
surface  current  require  high  harmonic  content  (p  large)  which,  from  (1-22), 
require  very  small  skin  depths  (high  conductivity)  for  the  high  frequency 
skin  effect  approximation  to  be  valid. 

In  Fig.  1-11  the  resistance  R p/^D  **  plotted  against  the  number  of 

harmonic  terms  used  in  the  series  representing  the  current.  The 

coefficients  obtained  by  either  of  the  approximate  methods,  unlike 

file  Fourier  coefficients,  are  a  function  of  the  number  of  terms  q  used 

in  the  series.  They  approach  the  exact  coefficients  in  the  limit  as  q 

becomes  large  or,  in  terms  of  the  resistance,  as  R^/R^  converges  to  the 

limiting  value.  For  this  reason  the  coefficients  used  in  constructing 

Fig.  1-11  are  those  found  for  the  limiting  value  of  R^/Rq.  From  Fig. 

1-11,  6  harmonic  terms  are  sufficient  to  give  R  /R_  to  within  1%  of  the 
>  p  o 

limiting  value  for  the  minimum  conductor  spacing s  presented  in  Table  1-1. 
Using  equation  (1-22)  with  6  harmonic  terms,  the  high  frequency  skin 
effect  approximation  will  be  valid  provided 


1 


2 


3  4  5 

Number  of  Harmonic  Term* 


6 


FIG.  1-11  THE  RESISTANCE  Rj/Ro  AS  A  FUNCTION  OF  THE  NUMBER  OF 
HARMONIC  TERMS  USEO  TO  REPRESENT  THE  CURRENT 
DISTRIBUTION  ON  EACH  WIRE.  THE  SPACINGS  C/o  ARE  THE 
MINIMUM  VALUES  PRESENTED  IN  TABLE  1-1. 
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a/d  »  1  (l-69a) 

8 

d  .* 

9(~)2 

-  g-  «  1  (l-69b) 

d-9^) 

These  conditions  are  satisfied  by  most  wire  sizes  used  in  practical 
antennas  operating  at  frequencies  above  1  MHg. ;  for  example:  l/8  inch 
radius  copper  wire  has  die  following  values 


Frequency  (MHg. ) 

a/dg 

(3dg/a)2 

(l-9dg/a) 

1 

34 

1.  1  x  10"2 

10 

107 

8.  5  x  10‘4 

100 

340 

8.  Ox  10"5 

The  failure  of  the  skin  effect  approximation  for  extremely  close 
conductor  spacings  places  no  serious  restriction  on  the  usefulness  of  the 
solution,  since  in  practical  applications  the  minimum  spacing,  deter¬ 
mined  by  the  thickness  of  the  wire  insulation,  is  usually  within  the  range 

of  values  covered  in  Fig.  1-9  and  Table  1-1  . 

% 

4.  Comparison  with  the  Work  of  Butterworth 
As  previously  mentioned,  Butterworth  has  calculated  the  ohmic 
resistance  of  systems  of  parallel  wires.  In  this  section  two  of  his 
formulas,  rewritten  in  the  form  R^/R^,  are  compared  with  the  present 

/ 
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calculations  of  additional  ohmic  resistance  due  to  the  proximity  effect 
(Fig.  1-9).  The  first  graph  in  Fig.  1-12  is  a  comparison  with  Butter- 
worth's  "semi-empirical  formula"  which,  for  high  frequencies,  can  be 
written  as  [10,  p.  709,  equation  53] 


R 

_JE  i 
R_ 


1 

8  wn<l/c) 


4 


un(a/c)' 


2(l-^vn(a/c)V 


(1-70) 


where  u  ,  v  ,  and  w  are  constants  which  depend  on  the  number  of 
n'  n'  n 

conductors  in  the  system.  This  formula  gives  results  which  are  in  fair 
agreement  with  the  present  calculations. 

In  the  second  graph  of  Fig.  1-12,  the  present  theory  is  compared 
with  another  of  Butterworth's  formulas,  one  which  is  often  found  in 
handbooks  on  coil  design  [11],  [13],  and  [14].  This  formula  is  derived 
by  making  assumptions  similar  to  those  already  discussed.  Consider 
each  conductor  to  be  in  a  uniform  magnetic  field  due  to  the  other 
conductors.  With  (1-53)  and  (1-4),  the  power  loss  per  unit  length  in  the 
m**1  conductor  is 


Watts/meter 

(1-71) 


and  the  resulting  ohmic  resistance  due  to  the  proximity  effect  becomes 


\  (a/c)2 


n 


2 
m=  1 


jf/m 


(1-72) 


FIG1-I2  THE  ADDITIONAL  OHMIC  RESISTANCE  DUE  TO  THE 
PROXIMITY  EFFECT-COMPARISON  WITH 
BUTTERWORTH'S  SOLUTIONS. 
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As  seen  in  Fig.  1-12,  this  formula  gives  results  which  are  obviously  in 
error  for  spacings  in  the  range  1  <  c/a  <  2.  It  is  applicable  only  in  a 
region  {c/a  »  1)  where  the  proximity  effect  is  of  little  interest,  since 
values  of  R^/Rq  are  small  and  fairly  independent  of  the  number  of  wires. 


5.  Optimum  Conductor  Spacing  when  the  Cross  Sectional 
Dimensions  are  Restricted 

In  certain  applications  a  given  number  n  of  parallel,  in-line 
conductors  must  fit  within  a  specified  length  i  ;  see  Fig.  1-13.  It  is  of 
interest  to  ask  for  which  wire  radius  a,  or  spacing  c/a,  is  the  resistance 
of  the  wires  a  minimum.  If  there  were  no  proximity  effect,  making  the 
radius  of  the  wire  as  large  as  possible  (a  =  t  /2  n)  would  minimize  the 
skin  effect  resistance.  With  the  proximity  effect  present,  increasing 
the  wire  radius  increases  the  loss  due  to  proximity  and  a  minimum 
resistance  point  is  reached  where  the  decrease  in  skin  effect  loss  is 
just  balanced  by  an  increase  in  proximity  loss.  In  Fig.  1-13,  the 
dimensionless  quantity  Zi r  i  R/nRS,  which  is  proportional  to  the  ohmic 
resistance  per  unit  length  of  the  system  of  conductors,  is  plotted  against 
the  normalized  wire  radius  a/f  .  The  points  of  minimum  resistance  are 
clearly  exhibited  in  Fig.  1-13  and  the  corresponding  conductor  spacings 
are  listed  in  Table  1-2. 


Number  of 
Conductors  n 

a/i 

c/a 

2rrfR 

Rs 

2 

0.  250 

1.00 

5.  33 

3 

0.  148 

1. 19 

10.41 

4 

0.098 

1.  37 

16.  07 

5 

0.  071 

1.  50 

22.01 

6 

0.  056 

1.  59 

28.  10 

n 

• 

0.  046 

1.  66 

34.  30 

8 

0.039 

1.71 

40.  57 

Table  1-2.  Conductor  Spacings  for  Minimum  Resistance. 


FIG.H3  THE  OHMIC  RESISTANCE  AS  A  FUNCTION  OF  THE 
WIRE  RADIUS  a  WITH  THE  DEPTH  OF  WINDING 
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6.  Conclusion 

Systems  of  equally  spaced,  in-line  conductors  carrying  equal 
currents  in  the  same  direction  have  been  studied.  A  set  of  integral 
equations  was  formulated  to  determine  the  transverse  distribution  of 
axial  current  at  high  frequencies  when  die  current  is  confined  to  a  dun 
skin  near  the  conductor  surface.  Using  the  integral  equations,  an 
approximate  solution  for  the  current  in  the  form  of  a  trigonometric 
series  was  obtained.  For  two  wires,  the  approximate  solution  for 
the  current  showed  good  agreement  with  an  exact  expression  obtained 
by  a  conformal  mapping  procedure. 

With  the  current  distribution  determined,  the  high  frequency 
resistance  per  unit  length  of  the  system  was  calculated  for  various 
numbers  of  conductors  and  spacings.  The  results  of  these  calculations 
may  be  summarized  qualitatively  as  follows: 

i.  For  small  numbers  of  conductors,  die  additional  ohmic 
resistance  due  to  the  proximity  effect  Rp/RQ  increases  either  with  an 
increase  in  the  number  of  conductors  or  with  a  decrease  in  the  conductor 
spacing.  This  was  checked  for  systems  with  up  to  eight  wires  and 
spacings  as  close  as  c/a  =  1.  1. 

ii.  For  closely  spaced  conductors  the  additional  ohmic  resistance 
due  to  the  proximity  effect  can  be  greater  than  the  resistance  of  the 
isolated  wires. 

iii.  When  the  cross  sectional  length  1  =  2  a  +  (n-1)  c  of  the  group 
of  conductors  is  restricted,  there  is  a  definite  wire  radius  that  will 
give  a  minimum  resistance  per  unit  length  for  the  system. 
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Only  cylinders  carrying  equal  currents  in  the  same  direction  were 
considered  in  this  chapter.  With  a  simple  scaling  of  the  harmonic  terms 
on  each  conductor  ii/e  present  theory  and  associated  computer  codes 
could  handle  systems  of  wires  with  different  net  currents  in  each  wire. 
Such  a  solution  would  be  useful  for  making  computations  for  multi  wire 
transmission  lines  where  the  wires  carry  currents  with  equal 
magnitude  but  in  opposite  directions. 


SECTION  n 


THE  ELECTRICALLY  SMALL  MULTITURN  LOOP  ANTENNA 

L  Introdaction 

The  single  tarn  loop  has  been  Hie  subject  of  much  investigation 
and  from  the  practical  standpoint  adequate  design  data  are  avail- 
able  for  this  structure  [28],  [29).  The  multiturn  loop,  with  no 
restrictions  on  electrical  size,  has  received  much  less  attention. 
The  solutions  available  are  for  the  "one  dimensional"  current 
distribution  and  therefore,  strictly  speaking,  only  valid  for 
loops  with  spacings  between  turns  large  compared  to  Hie  wire 
diameter  [30],  [31 J. 

In  practical  applications,  the  electrically  small  loop  ie 
often  used  because  it  has  a  desirable  field  pattern  as  compared 
to  larger  loops  whose  patterns  have  many  lobes.  The  ohmic 
resistance  of  small  loops  is  in  general  much  larger  than  the 
radiation  resistance,  thus  radiation  efficiencies  are  very  low 
and  greatly  dependent  on  the  ohmic  resistance.  In  an  effort  to  in  - 
crease  the  radiation  efficiency  multiturn  structures  are  often 
used.  The  radiating  properties  (radiation  resistance  and  field 
pattern)  of  electrically  small  single  or  multiturn  loops  are 
easily  derived,  either  directly  from  the  integral  form  of  Maxwell's 
equations  [28],  [32]  or  as  a  limiting  case  of  one  on  the  more 
general  analyses  mentioned  above.  These  methods  are  usually 
concerned  with  perfectly  conducting  wires  and  thus  provide  no 
information  about  ohmic  loss  of  the  antenna. 
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The  ohmic  resistance  of  a  small  loop  is  usually  taken  to 
be  the  same  as  mat  of  an  equivalent  length  of  straight  conductor. 

This  assumption,  although  adequate  for  the  single  turn  loop,  is 
not  for  the  multiturn  case.  In  a  multiturn  loop,  the  distribution 
of  current  over  the  conductor  cross  section  is  determined  by  the 
same  effects  discussed  in  Chapter  I  -  -  proximity  and  skin  effects. 
The  increase  in  ohmic  resistance  due  to  the  proximity  effect, 
which  is  normally  unimportant  in  large  antennas,  has  a  dramatic 
effect  on  calculations  of  the  power  radiated  by  electrically  small 
transmitting  loops. 

2.  Review  of  Small  Loop  Theory 

The  properties  of  electrically  small  loop  antennas  covered 
in  the  literature  are  briefly  discussed  below.  For  a  more 
detailed  discussion,  see  King  [32]  or  King  and  Harrison  [28]. 

The  model  chosen  to  represent  the  multiturn  loop  antenna 
is  illustrated  in  Fig.  2-1.  All  turns  of  the  loop  are  circular  and 
tie  in  parallel  planes.  The  straight  segments  of  wire  interconnecting 
the  turns  and  the  feed  wires  of  the  delta- function  generator  are 
short,  parallel  and  closely  spaced.  These  are  assumed  to  have 
negligible  ohmic  resistance  compared  to  that  of  the  overall 
circuit,  and  to  contribute  negligibly  to  the  radiation  resistance 
since  parallel  segments  carry  equal  and  oppositely  directed  currents. 
The  dimension  2c  is  exaggerated  in  Fig.  2-1. 

The  multiturn  loop  with  n  turns  will  have  essentially  the 
same  total  current  (I)  at  any  conductor  cross  section,  provided 
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the  total  length  of  the  loop  is  much  less  than  the  free  space  wave¬ 
length  at  the  operating  frequency.  More  specifically. 


I(s)  =1 


when 


(2-1) 


0onb<<l 

For  this  analysis,  the  following  additional  constraints  are  placed 
on  the  wire  radius  a  and  the  turn  spacing  c. 


a«b  ,  Poa<r<l 


2  2  2 
n  c  «b 


c  a 


A  real  power  equation  expressed  in  terms  of  the  scalar  and 
vector  potentials  $  and  A  for  the  loop  antenna  is 

.2 


J-  E  dv 

O 


)  * Re  ( 1 4- dv  - iu  j 


dv  -  iw  /  J*  A  dv  -  iu 

v  s 


/  H  Ids) 


(2-2) 


(2-3) 

(2-4) 


(2-5) 


where  J  represents  the  free  current  density,  n  the  free  surface 

charge  density,  a  the  conductivity  of  the  antenna  wire,  and 

the  electric  field  of  the  delta- function  generator.  The  first  three 

integrals  are  over  the  volume  occvipied  by  the  loop  conductor  and 

generator  while  the  fourth  is  over  the  surface  of  this  volume.  An 

e  1W*  time  dependence  is  assumed. 

Fig.  2-2  shows  sections  of  two  typical  loop  turns  and  the 

coordinates  associated  with  them.  Making  use  of  (2-1),  and 

assuming  the  transverse  current  distribution  to  be  the  same  at 

til 

any  cross  section,  the  current  density  on  the  m  turn  becomes 


FIG.  2-2  SECTIONS  OF  TWO  TYPICAL  LOOP  TURNS 
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Jm(r,  9,  q>)  =  Jnw(>'.9><P 


where 


tt  a 

/  / 

*9  =-rt  r=0 


(r,  9)  rdrdQ  =  I 


The  assumption  of  constant  current  also  precludes  the  possibility 
of  a  charge  accumulation  on  the  loop  turns.  Therefore, 


n  =o 

With  (2-6),  (2-7)  and  (2-8),  and  the  definition  of  the  delta- function 
generator 


=  -  V06(s) 

(the  distance  s  is  shown  in  Fig.  2-1),  (2-5)  becomes 

tt  „  a 

n  [ 

V  1  ?.TT  V. 

IV 


/  tt  „  a 

lH  / 


2 

J  (r,  o  )  rdrdo 
mcp 


a=-n  r=0 
tt„  a 


+  Re  [  -i2nbw 


'mtp  <r-  “>  Anfr  9>  rdrd9l 


9=-tt  r=0 


} 


Am^  Is  the  component  of  the  vector  potential  tangent  to  the  axis  of 

tVl 

the  conductor  of  the  m  turn. 

Referring  to  Fig.  2-2,  the  vector  potential  at  point  A  due 


(2-6) 


(2-7) 


(2-81 


(2-9) 


(2-10) 


to  the  current  element  at  point  B  is 


-53- 


iS  H 

e  ’■'o  mjt 

dA  (r,  a,  cp)  =  -x-k -  j£  (r',  0')  (b  +  r'  cos  ) 

micp  mi.  ^ 


cos  (Q  -  S'  )  r'di?  d0'dcp' 


(2-11) 


where 


Rmjl  =  f4b2®in2  [  (0  “  9*  )/2  ]  +  4(m-42  c2  +  r2  +  r'2 
+  2rr'  cos  -  cp'  )^  1 

The  vector  potential  at  A  due  to  the  current  in  ail  turns  is  then 

TT  TT  a 


A 


mcp 


(r‘,  0') 


eifoRmA 

Rmjfc 


cos  (0  -  q')  (b  +  r'  coscp')  r'dr' 


(2-12) 


(2-13) 


Introducing  the  condition  on  the  Conductor  length  described  in 
(2-2),  the  exponential  in  the  integrand  of  (2-13)  is  expanded  in 
a  power  series  in  0o^mJp  Keeping  the  first  two  imaginary  terms 
in  this  series  yields 


Im 


(eie°Rm*)  = 


Ko  mi, 


3  3 

3  R  . 

Ko  mi 


(2-14) 


and  the  imaginary  part  of  (2-13)  becomes 

n  TT_  a 


Ini(A  ) 
mcp 


’*■  I1/  I 

cp  s-tt  q  s-tt  r'=0 


3  _  2 
B  R  p 
o  mi 

- g - 


cos(0-  )(b  +  r*  cos  cp')  r'dr'd0  'dcp' 


(2-15) 
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For  the  purpose  of  calculating  the  second  integral  in  (2-10),  the 

radiation  term,  and  approximate  value  of  the  vector  potential 

A_  is  used.  Subiect  to  the  restrictions  on  the  conductor  radius 
mep 

expressed  in  (2-3),  is  approximately  the  vector  potential 

that  would  exist  on  the  surface  of  the  wire  with  the  loop  current 
I  located  along  the  axis  of  the  conductors.  King  [33]  discusses 
the  validity  of  this  type  of  approximation  when  used  in  calculating 
the  vector  potential.  With  this  simplification,  (2-15)  reduces 
to  the  following 


Im(A  ) 

mep 


u  rb 

4n 


cos(cp -  cp')dcp' 


(2-16) 


where 


R 


=  j4b2  sis 


ml  / - in<*  f  <0-  fl')/2  ]  +  4(m-£)2  c2  +  a2 

and  (2-10)  becomes 


j 


2  fib 

a 


J  (r,  a  ) 
mg'  ’  q  ' 


rdrdfl 


J  V  fy. 

a=-n  r=0 


fl3  l2 
l*  9  b 
,  o  o 

h — it - 


n 

y 

7=1 


TT 


V=-n 


-  RmjJ  cos(o  -  q  1  )d©  * 


(2-17) 


(2-18) 


Evaluating  the  second  integral,  (2-18)  becomes 
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vo  ■  ifi  [*■  f  f  — .]  ♦  **■*■*, y> 
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=  1  rOhmic  +  RRad. 


where  the  two  terms  on  the  right  hand  side  of  the  equation  are 
identified  as  the  ohmic  and  radiation  resistances. 

The  radiation  efficiency  of  the  n  turn  loop  is  now 


E  =  Rad. 

A  R„  ,  +  Rn,  . 

Rad.  Ohmic 


2  2  4,  4 
20rr  n  0Qb 


+Rohmic 


This  simple  form  is  a  consequence  of  the  constant  current 
assumption  which  makes  the  ohmic  and  radiation  resistances 
appear  as  circuit  elements  in  series. 

3.  Transverse  Current  Distributions 

To  evaluate  the  expression  (2-20)  for  the  radiation 
efficiency  the  transverse  current  distribution  is  needed.  If 
the  skin  effect  approximation  applies,  a/ds  >>  l,  the  ohmic 


resistance  term  in  equation  (2-20)  can  be  replaced  by 
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where  )  i*  normalized  surface  current  density  on 

equivalent  perfectly  conducting  loop. 

Using  a  procedure  similar  to  that  in  section  1-3,  the 

transverse  current  distribution  g  (ft)  can  bn  derived.  The 

emcp'0' 

integral  for  the  vector  potential  component  A  at  a  point 

th  * 

(r,  0,  cp)  just  off  the  surface  of  the  m  turn  is 
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where 


Rmf  2l>2  +  c2  +  r2  +  a2  +  4(m-£)c  (r  cosq  -  a  cosg  ') 

-  2b  (r  sing  +  a  sing  ')  -  2ar  cosgcosg'J  -  2[  (b  -  r  sing)  (b  -  a  sing  ) 
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compared  to  unity.  Dropping  terms  of  this  order,  (2-22)  becomes 


cos(tp-  q$) 

■  . 

q  -  p  cos(cp-  eff)  2 


dg'dcp* 


(2-24) 


This  is  equivalent  to  considering  the  quasistatic  fields  as  the 
primary  factor  in  determining  the  transverse  current  distribution. 
The  integration  with  respect  to  cd  may  be  expx cased  in  the  form 
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where  K  and  E  are  the  complete  elliptic  integrals  of  the  firs^  and 
second  kind  [34],  The  modulus  k  and  complimentary  modulus  k* 
of  the  elliptic  *ntegral*  £re 


=  QZ? 

sp  +  q 


(k')2  =  1  -k2 


Subject  so  the  restrictions  imposed  S/(Z-Z)  and  (2-4) 
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Since  rmjJ  is  of  the  order  of  2(m-l)c,  (k1)  is  a  small  quantity 


V)2 


(k*)  50|^|  «1 

and  the  power  series  representations  for  K  and  E  are  useful  [35], 
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Subst  eating  ine  above  serifs  (2-25)  ana  dropping  all  terms 

2  2  2 

sr- ^11  cornr>ared  to  unUy  &/b,  (m-l)  c  /b  or  less,  the  integral 
r  jr  the  vector  pot',  ati-'j  <  becomes 


Except  for  a  term  with  only  z  dependence  and  an  additive  constant, 
this  expression  is  the  same  as  that  for  A  (r,  o,  z)  in  the 
equivalent  system  of  parallel,  straight  conductors,  equation  (2-30). 

Due  to  the  symmetry  already  assumed  in  this  problem, 
only  /he  A  component  of  the  vectoi  potential  is  involved  in  the 
Py mss/'!  the  surface  current  density,  P.n(5  h 


«In(6) 


2AnMjr’  9  •  ^ 


r=a 


With  (2-31)  substituted  into  (2-32),  the  resulting  equation  for  the 
current  density  is  identical  to  that  for  the  straight  conductors  (1-39). 

Subject  to  the  inequalities  presented  in  equations  (2-2),  (2-3) 
and  (2-4),  the  transverse  current  distributions  on  the  loop  turns 


(2-31) 


(2-32) 
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and  the  ohmic  resistance  per  unit  length  are  thesame  as  those  for 
a  system  of  parallel,  straight  conductors  which  have  the  same  wire 
radius  and  spacing. 


4.  Radiation  Efficiency 

With  the  results  of  the  last  section  and  equation  (2-20),  the 

radiation  efficiency  of  an  n  turn  electrically  rmall  loop  is 
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Rearranging  terms,  the  efficiency  becomes 


8.  48  x  10"10  yff 
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(2-34) 


where  a'  and  b'  are  the  radius  of  the  wire  and  the  radius  of  the 
loop  normalized  to  the  free  space  wavelength,  fj^j^  is  the 
frequency  in  megahertz,  and  ar  i  s  the  ratio  of  the  conductivity 

7 

of  the  loop  wire  to  that  of  copper  (acu  58  5.  8  x  10  /  ohm-m).  In 

Fig.  2-3,  the  efficiency  is  plotted  as  a  function  of  the  dimensionaless 
3  _ _ — . 

quantity  (b1)  a'  /  \| '  an<*  '^e  num^er  turns.  The  dashed 
tines  are  for  no  loss  due  to  proximity  (R^/Rq  =  0)  while  the  solid 
lines  include  the  proximity  effect  for  a  spacing  c/a  =  1.  10.  For 
most  practical  applications,  these  two  lines  will  give  an  upper 
and  lower  bound  on  the  efficiency  obtainable  with  various  turn 
spacings. 


Neglecting  the  proximity  effect  can  lead  to  large  errors  in 
the  calculation  of  radiation  efficiency.  For  example,  from  Fig, 


Fft  E-J  RADIATION  EFFICIENCY  OF  ELECTRICALLY  SHALL  LOOPS 
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2-3,  without  the  proximity  effect,  the  calculated  efficiency  of 
a  three  turn  loop  can  be  larger  than  the  actual  efficiency  of  an 
eight  turn  loop  of  the  same  size  with  close  conductor  spacing 
(c/a  =  1.  10).  When  the  loop  is  used  as  a  transmitting  antenna 
the  radiated  power  is  directly  proportional  to  the  radiation 
efficiency.  Neglecting  the  proximity  effect  can  make  the  computed 
efficiency  for  a  small  loop  in  error  by  a  factor  of  two  or  larger; 
thus  errors  in  the  calculation  of  radiated  power  can  be  as  large 
as  one  hundred  percent. 

In  some  applications  a  constraint  is  placed  on  the  volume 
the  loop  antenna  can  occupy.  If  the  depth  of  winding  l  is 
restricted  to  a  value  much  smaller  than  the  diameter  of  the 
loop  (j i  «  b,  see  Fig.  2-4)  the  results  of  section  1-5  can  be 
used  to  optimize  the  efficiency.  With  no  proximity  effect,  the 
maximum  efficiency  is  obtained  when  a  =  £/2n  and  is  independent 
of  the  number  of  turns  n. 
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(2-35) 


If  the  proximity  effect  is  included,  the  antenna  has  optimum 
efficiency  when  the  turn  spacings  ire  those  presented  in  Table 
1-2 

E  =  _ _ i _ i _ 

A  ,  +  8'48lW'l0Jg 

n(b')5  J'  '  Ps  / 
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where  i'  =  f/V.  and  (2tifR/RS)  is  the  value  given  in  Table  1-2. 

Both  equations  (2-35)  and  (2-36)  are  graphed  in  Fig.  2-4, 

As  the  number  of  turns  is  increased,  the  term 

in  (2-36)  increases,  causing  a  decrease  in  the  efficiency. 

With  the  antenna  restricted  to  a  volume  of  this  shape,  it  is 
better,  then,  to  optimise  the  wire  size  rather  than  to  increase 
the  number  of  loop  turns.  At  a  first  glance,  this  last  statement 
seems  contrary  to  the  common  notion  that  increasing  the  number 
of  loop  turns  increases  the  radiation  efficiency.  It  must  be 
kept  in  mind  that  one  usually  speaks  of  increasing  the  number  of 
turns  while  keeping  the  wire  radius  and  spacing  constant,  some¬ 
thing  which  is  impossible  to  do  when  the  depth  of  winding  f  is 
also  fixed. 

Power  is  usually  supplied  to  electrically  small  antennas 
through  a  suitable  matching  network.  The  components  in  the 
matching  network  often  introduce  losses  as  large  as  the  ohmic 
loss  of  the  antenna.  The  overall  radiation  efficiency  of  the 
antenna- matching  network  combination  is  then 

E  *  EA  *  EM  •  <2*37> 

where  and  are  the  efficiency  of  the  antenna  and  matching 
network  individually.  In  this  chapter,  only  E^  is  considered; 
for  a  discussion  of  matching  network  efficiency,  see  Wheeler  [36], 


5. 


Conclusion 


The  analysis  in  this  chapter  has  shown  that  the  results  obtained 

for  the  ohmic  resistance  per  unit  length  of  a  system  of  straight 

wires  are  applicable  to  the  electrically  small  multiturn  loop  when 

the  depth  of  winding  of  the  loop  is  small  compared  to  the  loop 
2  2 

radius,  (nc)  <  <  b  .  Two  separate  calculations  of  the  radiation 
efficiency  of  small  multiturn  loops  were  made:  the  first  includes 
the  added  resistance  due  to  the  close  proximity  of  turns  and  the 
second  neglects  all  proximity  losses,  i.  e.  considers  the  ohmic 
resistance  of  the  loop  to  be  the  same  as  that  for  an  equivalent 
length  of  straight  conductor.  A  comparison  of  the  results  for 
these  two  cases  indicates  that  the  proximity  effect  is  an  important 
factor  in  making  accurate  calculations  of  radiation  efficiency, 
especially  for  loops  whose  efficiency  is  below  10%. 

The  problem  of  optimizing  the  radiation  efficiency  of  an 
electrically  small  loop  confined  to  a  fixed  volume  was  also 
examined.  The  special  case  of  a  circular,  multiturn  loop 
restricted  to  a  volume  whose  depth  is  small  compared  to  the 
loop  radius  {£  c<b)  was  treated.  For  this  geometry  there  is 
an  optimum  wire  radius  which  gives  maximum  radiation  efficiency 
for  a  given  number  of  loop  turns.  With  the  optimum  wire  radius 
used  for  each  number  of  turns,  the  radiation  efficiency  was 
found  to  decrease  with  an  increase  in  the  number  of  turns,  indi¬ 
cating  that,  from  the  efficiency  standpoint,  it  is  better  to  optimize 
the  wire  radius  than  to  increase  the  number  of  turns. 


The  change  in  the  transverse  distribution  of  current  due 
to  the  proximity  effect  wiij  also  alter  the  loop  inductance.  The 
inductance  of  the  loop,  however;  does  not  have  to  be  known  to 
a  high  degree  of  accuracy  in  most  applications,  since  it  is 
usually  made  to  resonate  with  a  variable  capacitance  in  a 
matching  network. 


SECTION  m 


EXPERIMENTAL  INVESTIGATION 

L  DeitripKoB  of  Btperiinent*!  Apparatus 

To  verify  the  results  of  Section  I.  experimental  apparatus 
was  constructed  for  measuring  the  transverse  distribution  of 
current  on  a  system  of  parallel  round  wires,  see  Figs.  3-1  acd  3-2. 
The  parallel  wires  are  modelled  fay  34“  long  copper  tubes  inter¬ 
connected  with  wire  braids  so  that  they  carry  equal  currents  in 
the  same  direction.  A  100  Watt,  100  KHx  transmitter  drives  a 
current  of  the  order  of  1-2  Amps,  through  the  model,  which  for 
matching  purposes  is  fed  in  series  with  a  US  Ohm  load.  The 
current  distribution  is  measured  by  sampling  the  transverse 
magnetic  field  with  a  small  loop  probe  mounted  on  one  of  the 
tubes.  This  tube  has  plugs  fitted  with  beryllium  copper  finger 
stock  at  both  ends;  these  maintain  electrical  contact  as  the  tube 
is  rotated  (Fig.  3-3  ).  The  voltage  at  the  terminals  of  the 
loop  probe  is  amplified  and  metered  using  a  General  Radio 
model  1232-A  Toned  Amplifier  and  Null  Detector. 

At  100  KHz  the  1^"  copper  pipes  are  about  200  skin  depths 
in  diameter;  thus  the  axial  cor  rents  are  confined  to  a  thin 
layer  near  the  outside  of  the  tube.  The  tubes  are  also  about 
20  skin  depths  thick,  so  they  are  electrically  equivalent  to  solid 
conductors. 

To  maximize  the  angular  resolution  of  the  measured 
current  distribution,  the  radial  dimension  of  the  loop  probe 
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FIG.  3-1  BLOCK  DIAGRAM  OF  THE  CURRENT  MEASURING  EQUIPMENT 


-71- 


was  made  as  small  as  possible  (0. 050").  Sli  ce  the  fields  are 
fairly  uniform  over  small  lengths  near  the  center  of  the  tube, 
fee  axial  dimension  of  the  loop  coold  be  a  few  inches  long.  Using 
die  theorv  of  Wbiteside  and  King  (37]  die  voltage  at  the  terminals 
of  die  rectangular  loop  when  the  tube  carries  a  total  ccrrent  of 
one  Ampere  is 
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where  I  ,  |  are  the  radial  and  axial  dimensions  of  die  locp  in 
*  « 

inches,  r  the  loop  wire  radios  in  inches,  and  the  load  impedence 
at  die  probe  terminals  which  is  about  50  K  Ohms  for  die  G.  R. 
U.S?«A.  From  (  3-1  t,  die  3”  x  0.  050"  loop  probe  constructed 
of  28  gage  wire  provides  a  0. 6  m  volt/Amp  signal.  This  is  more 
than  atf e&itn ^  far  metering  on  die  G.  R.  1232-A,  since  it  hat  a 
maximum  sensitivity  of  10  u  Volts  for  a  full  scale  deflection  at 
100  KHz,  For  rigidity  a  polyfoam  support  was  placed  between 
die  loop  and  the  tube  (see  Fig.  3-3  >. 

Initial  measurements  indicated  that  the  metering  circuit 
wan  picking  up  a  very  strong  signal  induced  by  the  large  loop 
formed  by  the  tub?-*  and  interconnecting  wires.  To  eliminate 
this  interference,  die  meter  was  completely  enclosed  in  a 
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copper  box  and  all  cables  used  were  doubly  shielded. 

The  G.  R.  1232-A  Toned  Amplifier  and  Null  Detector  was 
calibrated  at  100  KHz  using  a  pair  of  Hewlett  Packard  precision 
attenuators  as  a  standard.  Fig,  3-4bis  a  schematic  of  the  circuit 
used  for  the  calibration.  The  linear  scale  meter  reading  is 
plotted  against  the  attenuator  setting  in  Fig.  3-4a.  The  small 
vertical  lines  indicate  the  experimentally  determined  points; 
a  +  unit  reading  error  is  assumed.  The  solid  line  was  constructed 
by  fitting  polynomials  to  the  experimental  points  over  three  ranges. 
The  polynomials  in  the  form  csed  to  correct  die  experimental 
data  are 


V*  sr  V  ,  63  <  V  <  100 
m  m  —  m  — 

.-2. 


V*  =V  +  3.  5-4.  8  x  10'^(V  -161-5. 7  x  10“4(V  -16)2,  16  <V  <63 
m  m  m  m  m 

V  =  V  4  3.5  .  0<  V  <16 
m  m  —  m  — 

where  V  is  the  meter  reading  and  V*  die  corrected  meter 
m  nn 

reading.  This  correction  is  apparently  only  a  function  of  the 
meter  circuitry  and  not  die  linear  amplifier  section  of  the 
instrument,  since  the  same  correction  applies  over  a  40dB.  range 
of  amplifier  gain. 

Since  the  theory  predicts  both  positive  and  negative  currents 
on  closely  spaced  tub  a  method  was  devised  to  experimentally 
verify  a  180°  phase  shift  in  the  current  density.  Referring  to 
Fig.  3-1  ,  a  small  terricc  loaded  loop  is  used  to  sample  the  field 
of  the  large  loop  formed  bv  the  tubes  and  the  interconnecting  wires. 
This  reference  signal  is  added  to  the  signal  from  the  current  probe 
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in  a  resistive  summing  network  and  then  metered.  The  phase  of 
the  current  on  the  tube  is  determined  by  noting  if  the  signal  from 
the  current  probe  adds  to  or  subtracts  from  the  reference  signal. 


2.  Correction  for  Interconnecting  Wires 

In  addition  to  the  net  currents  in  the  tubes,  three  other 
current  elements  influence  the  current  distribution  on  the  tube 
cross  section.  They  are  currents  in  the  horizontal  and  vertical 
interconnecting  wires  and  negative  line  currents  which  represent 
the  absence  of  a  continuation  of  the  axial  current  beyond  the  ends 
of  the  tube.  Referring  to  Fig.  3-5.  these  currents  can  be  treated 
as  filamentary  elements  since  each  is  at  a  distance  from  the 
probe  which  is  large  compared  to  the  tube  radius  (s  =  20a. 
w  =  60a).  Asa  result,  their  effect  on  the  transverse  distri¬ 
bution  of  current  is  additive  in  the  sense  that  it  may  be  subtracted 
from  the  measured  data  to  obtain  results  for  direct  comparison 
with  the  theoretical  distributions. 


The  vector  potential  at  a  point  (r,  0  ,  z)  near  the  center  of 


the  m*k  tube  is 
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Rmil  =  Rw“r*in6)2  +(2(m-^c  +  r  cosg)2  +  (z-a!  )2]2  (3-5) 

Rmi2  =  ^z'*z^2  +{2(m-i)c  +  r  cose)2  +(r  sin®)2  ]2  (3-6) 

Rmi3  =  ^*‘z^2  +<z<m-4c  +  r  cosa)2  +  (y-r  sing)2  ]2  (3-7) 

Rmj|l  s  f(s+z)2  +  (2(m-jQc  +  r  cosg)2  +  (y-r  sing)2  ]2  (3-8) 


is  the  normalized  surface  current  density  induced  in  the 
cylinders  by  the  three  external  current  elements.  It  is  the  term 
which  must  be  subtracted  from  the  measured  current  for  comparison 
with  theory. 

The  following  boundary  condition  relates  Am  and  gmC. 
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Substituting  (3-4  into  (3-S)  yields 
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where 


2  2  ~~ 

r<  =  [  (l-(a/w)sine)2 +(2(m-A)(c/w)  +(a/w)cosQ)  +(z‘/w)  ]2 

2  2  £  i 

r'  =  [(z'/a)2  +(2(m-jt)(c/s)  +(a/s)cosQ)  +(a/s)  (aine)  ]2 
mjt2 

r*  =  [  1  +  (2(m-jt)(c/a)  +  (a/ajcoag)2  +  (y/»  -  (a/a)ainQ)2  l2 
mi  3 

The  fir8t  integral  in  (3-10)  was  evaluated  in  8ection  1-3;  the  other 
integrate  are  a  atandard  form  f  1 9»  p.  50,  200.  03],  Performing 
theae  integrationa,  (3-10)  becomea 


g_,e(0 

me 


KmJe.e')  g^faW  -  2(a/w)(a/w) 


(3-11) 

(3-12) 

(3-13) 


(a/w  +  2(m-i)fc/w)  coaq  -  ainfl) _ _ _ ^ - 

[  (l-(a/w)ainflj  +  (2(m-A)(c/w)  +  (a/w)coao)£]  [  (i-(a/w)ainQ)  +  (2(m-$ 

,  .  2(a/s)  S  A^jLtBlSzmsI^L 

(c/w)  +  (»/w)coso)2  +(«/w)2]2  'z,  [  (2(m-f)(c/s)  +  (a/s) 

(.  I  _ 

coao)2+(a/a)2(8in9)  ]  [  1  +  (2(m-  ».)(«/»)  +  (a/ajeoag)  +(a/a) 

- r-ri  +2(a/a)ainq  V  — —  ,  - - - 2H; 

(ainfi)  ]2  j  ^  [  l  +  (2(m-jt)(c/a)  (a/a)coaf)  ] 

,  (w/a)(l  -  (a/w)  ainS _ _ _ _ + 

[  1  +  (w/a)2  (l-(a/w)aine)2  +  (2(m-jt)(c/a)  +  (a/a)coao)  ]2 
V 

1*1*)*M - 2 - j - 2-t-  I  (3-i4) 

[  1  +  (2(m-A)(c/8)  +  (a/8)co89)  +(a/8)  (8ino)  )2  \ 


When  terms  of  order  (a/s)  or  less  are  dropped,  equation  (3-14) 
reduces  to 


me 


•4-J  £ 

L.  .. 


jfm 


Kmf /«.*?•)  ♦*!»/•)  y 

i=l 


f  1  -fr(s/w)2  +  4(m-Jt)Z(c/,ar)^  ]*' 


(s/w)2fsino  -  2(m-l)(c/w)cosa  1 
[  I  +  4(m-4  (c/w)  ] 


sum 


[  1  +  4(m-jt)  (c/s)  J 


-  2{a/s)  y 

"  [(a/s)  +  4(m- f)  (c/s) 


+  4(m-  jO(c/s)(a/s)  J 


1  - 


[  i  +  4(m-j8Z(c/s)2  ]* 


(3- 


If  the  interaction  between  the  induced  currents  on  the  tubes 
is  ignored,  that  is,  each  tube  is  considered  as  isolated  from  the 
others,  the  integral  in  (3-15)  disappears,  and  a  first  order 
approximation  for  the  current  results. 


2(a/s) 


2 

(s/w)  fsinfl  -  2{m-j0(c/w)cosq  ] 

r  . .  a  m  it-.  i  i-  i 


[  1  +  4(m- 1)  (c/w)  ] 


sin9 

[  1  +  4(m-j&)^(c/8)^  ] 


-  2(a/s) 


(a/s  -f  2(m-j)(c/s)  cos| 

[(a/s)  +4(m-jJ  (c/s)  +  4(m-j9(c/s)(a/s)cos?  ] 
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I- 


'i* 


j  [  1  +  4fm-I|  (c/s)  ] 


€3-561 


The  three  terms  in  (3-16}  are  doe  to  the  current*  in  the  vertical 
interconnecting  wires,  horizontal  interconnecting  wires  and  axial 
tube  extensions,  respectively.  The  current*  induced  by  die 
horizontal  wires  are  the  major  factors  since  they  are  about 
(w/s)^  =  10  times  greater  than  those  due  to  the  vertical  wires  an’J, 
for  large  spacings  (c/s)  =  1,  at.  least  3  times  greater  than  those 
due  to  the  tube  extensions. 

To  solve  for  the  current  in  die  complete  expression  (3-15), 

which  includes  interactions  between  tubes,  a  trigonometric 

series  is  posutlated  for  g  (6).  Since  die  average  value  of 

g  (0)  is  zero  and  it  is  not  symmetric  about  die  lines  0:nor 
me 

tt/2  ;  3n/2,  the  series  has  no  constant  term  and  contains  both 
sine  and  cosine  terms. 


g  (6) 

8mc'  ' 


H 

5 


I  ‘mep  cos<^>  +  bmcp  sin<^>  1 


(3-17) 


Substituting  (3-17)  into  (3-15),  the  following  results  are  obtained 


n  even 


L 

P=i 


a  [  -  cos(pd)  +  • 

mep  1  '  tt 


n 


K{d  ,  d1)  cos  (pb')du'  ] 
m,  n+l-m 


H 

*5 


b  [  -  sin(pO) 
mep  1  'r 


.  Ml 


TT 


tt 

1 


K(6,  o')  sin(po')d0'  ] 
m,  n+l-m 


a  =-TT 
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n/2  q 

♦  E  E 

1=1  p=i 

l#m 


n/2  q  /• 

+  E  ^  blcp  4  J  [KJ^ »')  -  1-1)** K  ft  -0  1 


*cp 


nr  i 

I  t*  (9.  3*>  +{-0PK  (9.  0*}  ]  cos{poJ)daf  J 
-i*=-n  m*£  m.a+!-i  | 


o 

=  -2(a  /*) 


4=1 


■  . r~  ~ 

[l+(s/w>  +4{m-j&  (c/w)*] 


(»/w)Zfsin  -  Z(m-J|(c/w)co33i 

MM— 1  — 


[  1  +  4 (m-4  (c/w)  ] 


_ «inq 

TTT^Wi 


+  2(a/s) 


U 

E 


(a/s  -f  2(m-A)(c/a)co80 

.2  .  '  77F 


T  (a/s)  +  4(m-4£(c/s)  +  4{m-4(c/s) 


(a/s)  cos0  J 


_ 1 

[l+4(m-i)Z  (c/s)2  ] 


(3-18) 


m  =  1,  2,  -  —  n/2 
n  odd 


/pJ,  m  f  (n+L)/2  ) 
cos  j  ? 

(2pQ,  m  =  (n+l)/2 J 


K(«,  6')  cos(py')dB '  ) 
m,  n+l-m 


E 

P=i 


b  [  -  sin9< 
incp  L  > 


/p9  ,  m  $  (n+1)  /?. 

(^2(p+l)6,  m  =  (n+l)/2 
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(n-I)/2 


-l^h(TO)  /  K(0,*|  *iaii£')<&  ]  +  Y  V*  a, 

n  J  m,  n-H-m  JJ  Z  fcp 

5=1 


r 

/*—  /  [K  (S , 

(’  I—  -■* 

Kl  k."' 

•  ft 


8')  +  c-l)p  K  (S.  S')  ]  co*(pa’)di 
m,  n+1-1 


*)  (n-IJ/2 

A*  i 

J  i-i 


i=i 

Ifn 


)  q 

O')  -  (-1)P  K  (6  .  J)  ]  siatp&'Jde1  (  +  V  a 

m,  n+l-£  ) 

7  P=1 


nil  __ 

T  cp 


J 

i 

$  S«TT 


Tf 


K  <V  i\  „  <=«(2P9'>«'1  +  t  bm+l 

m,(n+x)/2  sl  “j — cp 


i¥*  J 

O’  =~TT 


K  (o  ,  Q*)  sin(2p+l)  s^dQ*  ]  =  -2(a/s)  Y 
m,(n+l)/2 


[  1  +  (s/w)  +4(m-8  (c/w)  ] 


(s/w)  fsjnfl  -  2(m-jQ(c/w)  coso]  +  sin8 

f  l  +  4(m-D  (c/w)  ]  [  1  +  4(m-jt)  (c/s)  ] 


n 


2(a/s)  y 


(a/s  +  2(m-l)  j^c/s)  C089) 


A=>1  [(a/s)  +  4(m-jQ  (c/s)  +  4(m-jQ(c/s)(a/s)  cosO] 


1 

1  -  . . . : — 5 - ’>~~T 

[  1  +  4(m-J t)  (c/s)  ]2 


m  =  1,  2,  -  —  (n+1) /2 


(3-19) 


The  definite  integrals  which  appear  in  (3-18)  and  (3-19)  can  be 
represented  by  two  forms.  The  integrals  containing  cos(p@?)  terms 
arc  the  same  as  the  integrals  denoted  by  J(c-  m-5.  pj  in  section 
1-5  and  evaluated  in  Appendix  A.  The  integrals  involving  siafpof) 
terms  are  of  the  form 

n 

I«(*.  ro-i.  p)  =  -L  /  1 ,1  *  -(c-,^  Co..' - 

15  /  [  4{m-l)  (c/a)  +  2  +  4(m-g)(c/a)cos» 

'a's-tt 

4  sinfl  licajj  SMx'pe^dg1 _ 

-  (4{m-i)(c.‘a)  4  2  cos(s'J)  cosS*  -  2  sing  sing1  ] 


p  =  1.  2 -  (3-20) 

An  evaluation  of  this  integral  is  in  Appendix  C,  the  results  tf  which 
are 

I’(S.  m-£.  p)  *  - -  [A*s2 +E's +C»  ]  (3-2la) 

(l-o  <-*)p 

where 


s 

2  2  - 
-  (4(m-l)  (c/a)  +1  +  4{m-i){c/a)  cosy)2 

(3-2!b) 

A* 

=  sin(a  -  (p-1)  t) 

(3-2lc) 

B' 

=  -2(1  +  2(m-l)(c/a)  cosa;  sin(p?) 

(3-21  d) 

C’ 

=  -sin(Q  +  (p+1)  *) 

(3-21e) 

an  *  (  -57 — ■  » ■■■*—  — -A  ,  (m- 1 )  -  1,  2,  - 

V  Z(m-2)(c/a)  tcosi)  / 

- 1  /  -sin°  \  ^ 

(  2(m-£)'(c/a)  +  cos0  j  .  (m- 1)  »  -l,  -2,  -  - 

The  principle  value  of  tan  ^  is  used  in  (3-2lf). 


Witt  {3-21),  the  system  of  cqaSoss  (3-18)  for  the  case  a  ma 


becomes 


q  n/2 

^  *mCJ>  [-  cos(pe)  +{-I)PIfS.  Zm-BB-1,  a)]  +  ^ 
p=l  £-» 


«5 


pel 


<9 

[Ife,  m-1,  ^  +  (-»P  IC#.  pil  *  J  [  -  s :*,*)-( -1) 

p=l 

n/2  q 

I*  (q  ,  2m-n-l,  p)]  +  ^  ^  bfcp  ^  **  ^ 

1=1  p=l 


-  (-i)Pl'  {5,  m+J-n-i,  p)]  =  -  2{a/s)  J 

IfI 


_  1 

[  1  +  (*/w)2  -5-  4{m-l)^(c/s»)^]^ 


(  f»/w)2f  sin9  -  ??m-l(c/w)  co sg 
(  [  l  +  4{m-j|  (c/w)  ] 
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S*mil»r  resales  are  obtained  for  the  case  n  odd.  Formula  (3-22) 
represents  a  set  of  a/2  eqoatioss  relating  the  qn  unknowns.  This 
system  of  equations  cac  be  solved  by  either  of  the  two  approximate 
methods  discussed  in  section  I>Sj,  methods  of  collocation  and  least 
squares. 

Appendix  D  contains  a  listing  of  a  computer  program  which 

solves  for  the  coefficients  a  and  b  by  the  method  of 

mep  mep  ' 

collocation.  The  qn  matching  points  in  the  internal  0^  6^  2^  are 
chosen  as 

Bk  * k  •  k  =  l- 2 "  ^ 

on  all  cylinders  except  the  center  cylinder  in  a  system  with  n  odd 
where  the  points  are 

9k  * 

Examples  of  the  correction  currents  g  are  plotted  in 

me 

Fig.  3-6.  Both  die  fall  correction  current  (3-17)  and  the  first 
order  correction  current  (3-15)  are  shown.  A  comparison  of  the 
two  curves  indicates  t’  -  '  the  interaction  between  induced  currents 
on  the  tubes  must  be  included  in  any  accurate  expression  fer  the 
correction  current. 

3.  Results  of  the  Measurements 

Current  distributions  were  measured  on  systems  with  up 
to  six  cylinders  and  spacing  s  ranging  from  cl  a  =1. 10  to  c/a  =  2.  50. 
After  correcting  for  meter  calibration,  the  measured  values  were 
normalized.  The  procedure  for  normalizing  was  first  to  measure 


(3-23) 


(3-24) 
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the  current  distribution  on  the  system  of  conductors  with  the  meter 
gain  held  constant  and  a  known  current  flowing  through  the  model. 
The  system  of  conductors  was  then  replaced  by  a  single  conductor 
and  the  current  distribution  measured  with  the  meter  gain  and 
current  through  the  model  the  same  as  in  the  previous  measurement. 
The  normalized  currents  on  the  multiwire  system  were  obtained 
by  dividing  the  measured  currents  by  the  average  value  of  the 
measured  currents  on  the  single  cylinder. 

The  measured  currents,  with  the  correction  current  g 

Iilv 

subtracted  cut  after  they  were  normalized,  are  compared  with 
the  theoretical  distributions  in  Figs.  3-7  to  3-11.  The  circles 
about  the  measured  points  indicate  the  range  of  error  (  +  2  scale 
units)  associated  with  the  repeatability  of  the  measurements.  The 
measurements  are  in  good  agreement  with  the  theory. 

The  minimum  apacings  used  in  making  the  measurements 
were  restricted  to  c/a  =  L  10  for  two  wires  and  c/a  *  1.  25  for 
three  or  more  wires.  For  three  or  more  wires,  the  currents  at 
adjacent  points  on  consecutive  cylinders  are  quite  large  when  the 
spacings  are  small.  The  radial  dimension  of  the  loop  probe 
is  a  significant  fraction  of  the  distance  between  cylinders;  for 
example,  when  c/a  =  1. 10  the  gap  between  the  cylinders  is  only 
three  times  as  large  as  £?.  Asa  result,  the  loop  responds  to  the 
currents  on  both  cylinders  giving  an  erroneous  interpretation  of 
the  current  density.  The  problem  is  not  as  severe  for  two  wires 
since  the  currents  at  adjacent  points  on  the  two  cylinders  approach 
zero  as  c/a  goes  to  l. 


FIG.  3-7  MEASURED  AND  THEORETICAL  SURFACE  CURRENT 
DISTRIBUTIONS  FOR  TWO  WIRES 


—  Theoretical 
O  Measured  points 


FIG.  3-8  MEASURED  AND  THEORETICAL  SURFACE 

CURRENT  DISTRIBUTIONS  FOR  THREE  WIRES 


T,T 


FIG.  3-9  MEASURED  AND  THEORETICAL  SURFACE 
CURRENT  DISTRIBUTIONS  FOR  FOUR  WIRES 


RG.  3-10  M5rA^dAN 


FIG.  3-11  MEASURED  AND  THEORETICAL  SURFACE 
CURRENT  DISTRIBUTIONS  FOR  SIX  WIRES 
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4.  Conclusion 

An  experimental  apparatus  was  constructed  to  measure  the 
transverse  distribution  of  current  on  systems  of  parallel  conductors. 
After  correcting  the  measured  data  for  equipment  calibration 
and  extraneous  sources,  good  correlation  between  theoretical 
and  experimental  current  distributions  was  obtained. 
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Using  the  trigonometric  identity  cos(A  E)  =  cos  A  cos  B  -  sin  A  sin  B 
to  combine  the  cos@'  and  sin@’  terms,  the  denominator  of  the  integrand 
becomes 


2 

D  =  r^^  ss  +  1  +  2s  cos(9'  +  ty) 

where 

2  2  - 
S  =  (4(m-£)  (c/a)  +  1  +  4(m-  A)(c/a)  cosq  )  2 


(A~2) 

(A -3a) 

(A -3b) 


The  principal  value  of  tan  *  is  used  in  (A -3b),  The  quantities  s  and 
are  identified  with  the  geometry  of  the  system  in  Fig.  A-l,  Applying 
standard  trigonometric  identities  the  numerator  of  the  integrand  is 
expanded,  giving 

N  =  I  A  cos[  (p-lj(?'  +  ty)]  +  B  cos  [  p(y  :  +  fy)] 

-  C  cos  [  (p+t)(o'  +([)]-  E  sin  [  (p- IXq’  +  \||)] 

+  F  sin  [  p(e 1  +  (i)]  -  G  sin  [  (p+l)(S  1  +  'll)]  j 


(A -4) 
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if 


I 


k 


where 


A 

=  cosfo  -  (p-1)  $  ) 

(A -5a) 

B 

=  2(1  +  3{m-i)(c/a)  coso)  cos(pf) 

(A -5b) 

C 

=  cos(0  +  (p+i)  $) 

(A-5c) 

E 

=  sin(g  -  (p-1)  $) 

(A-5d) 

F 

=  2(1  +  2(m-/)(c/a)  cose)  sin(p^) 

(A-5e) 

G 

=  sin(8  +  (p+l)$) 

(A-5f) 

In  terms  of  the  new  variable  £  =  equation  (A-l)  is 


I  (ft  ,  m-l, 


TT 

.  1  f  -A  cosf  (p  -  1)  4 1  + 

p>  *  -zr  |  - - 

I  s  *  1  +  2s  cos 

T  ="n 


B  cos  fp£l 

(4) 


-  C  cos  r  (p  -  l)4l  +  E  sin  fp£i  -  F  sin  f  (p  4  l)4ldn» 

The  8 in  [()$.]  terms  integrate  to  zero  and  the  remaining  terms  are 
in  the  form  of  a  definite  integral  which  is  readily  evaluated 
[  19,  p.  219,  858.  536] 


2tt 


C08 


LuSL 


1  +  H  cos 


=0 


ZTrfrjl-  l?  -1  ]p 

TO  L 

2ti 

,1  -  hV,J" 


,  P  =  1.  2 

,  p  =0 


Substituting  (A -7)  in  (A-6)  and  rearranging,  I(e,  m-l,  p)  becomes 

1 


I  (q  ,  m-l,  p)  = 


gjP"^  [AS  +B8+C],  P  =  1,  2  -  -  - 


-1 


8(1  -  S  ) 


[Rs  +  2C]  , 


p  =  0 


(A-6) 


(A -7) 


(A -8) 
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Appendix  B 

Listings  of  Computer  Programs 

This  appendix  contains  two  computer  programs  written  in 

Fortran  IV  language  for  use  on  the  I.  B.  M.  360/65  computer. 

Both  programs  compute  the  coefficients  of  the  trigonometric 

series  for  the  normalized  surface  current  densities,  g  ,  and 

m 

the  normalized  resistance  per  unit  length,  R/Rq.  The  input/ 
output  formats  for  both  programs  are  identical  and  specify 
the  following  parameters:  the  number  of  conductors  n,  the  number 
of  harmonics  q,  and  the  spacing  c/a. 


Preceding  page  blank 


s 


C  FORTRAN  IV  PROGRAM  FOR  SOLUTION  BY  THE  METHOD  OF  COLLOCATION 

C 

C 

C  THIS  PROGRAM  USES  THE  METHOD  UF  COLLCCATION  TO  SOLVE  A  SYSTEM  OF  EQUATIONS 
C  FOR  THE  COEFFICIENTS  OF  TRIGONOMETRIC  SERIES.  THE  SERIES  REPRESENT  THE 
C  NORMALIZED  SURFACE  CURRENT  OENSITY  ON  EACH  MIRE  IN  A  SYSTEM  OF  NW  EQUALLY 
C  SPACED,  PARALLEL,  PERFECTLY  CONDUCTING  HIRES.  USING  THIS  CURRENT  AN 
C  APPROXIMATE  VALUE  CF  THE  NOHPAlI ZED  UGH  FREQUENCY  RESISTANCE  OF  THE  SYSTEM 
C  IS  CALCULATED.  THE  NUMGER  CF  HARMONIC  TERMS  USED  TO  OESCRIBE  THE  CURRENT  ON 
C  EACH  WIRE  IS  NH.  THE  RATIO  CA  IS  EQUIVALENT  TO  THE  SPACING  BETWEEN  WIRE 
C  CENTERS  DIVIOEC  BY  THE  HIRE  DIAMETER.  FOR  NW  WIRES  AND  AND  NH  HAKMUNICS  THE 
C  SIZE  OF  TIE  MATRIX  TII.J)  MUST  BE  AT  LEAST  AS  LARGE  AS  T< NW1 *NH* l ,NW1 *NH) 

C  wHcRE  Nwl=NW/2  FOR  NH  EVEN  AND  NWl=SNW*l)/2  FOR  NW  CCD. 

LOGICAL  LSCLVE 

DOUBLE  PRECISION  P I , THETA ,sEP ,CA , SOURCE. T , SHI N 
CCMMCN  THETA, SEP, PI/MATR1X/T  (48, 4<J) 

WRl TEI6,l ) 

1  FORMAT  I IHI» 

PI*i.l4l5S2fc535897-»3D0 

2  READ ( 5,31  NW.NH.CA 

3  FORMAT (I  I, 3X, 12, 3X, Ob. 3) 

WR1TEI6.4)  NW.CA.NH 

4  FORMATI///35XIl,21H  WIRES,  SPACING  CA*,F6.i,4M,  ,12,10H  HARMON 

1ICS) 

NW20=12*Nw+U(-l)**(NW  +  l)  )/4 
Nw2E*I 2*Nw-lM-l)**NW)/4 
NSI ZE=NW2C*NH 
NAUG*NH*NW20+ 1 
DU  14  L 3 1 , NH 
NP.1  =  INW20-1)*NH+L 
C  SETTING  COLLGCATluN  POINTS 

THETA=Pl*IDFLOAT(L)l/(OFLOAT(NF+l)) 

I F( 2*1 NH/2) .NE.NH.ANU.L.LE. (NH*l>/2)  THETA* PI *FLO AT ( L) /FLOAT ( NH*2 ) 

I  FI  2*1  NH/2)  .NE.NH.ANO.L.GT.INH+U/2)  THETA3  PI  *FLUA  T{ L* l)  /FLUATC  NH* 

12) 

DO  14  P*L  «NR  1  «NH 

NKH  *  1 ♦  I M-L ) /NH 

IF  I Nw2Q. EQ.NW2E )  SO  TO  5 

IF  1NRW.EQ.NH20)  THETA* THETA/ 2. 000 

b  NSW  *  l 
N1=NH 

QG  14  N-l.NAUG 

IF  (N. EQ. NAUG)  CO  TO  11 

IF  IN-N1)  7,7,6 

6  N l* NI  ♦  NH 
NSw*NSW*l 

7  NSH=N-NH* (NSW-i  ) 

IF  INRW.EQ.NSWt  GO  TO  V 
SEP*CA*UFLOA  T (NSW-NRW ) 

IF  ( I  2  *NSW) .EQ.  INH* 1) )  GG  TU  8 
T (M,,N)*-SNININSH) 

SEP  *CA*OF  L.1AT  <NW*1  -NSW-NRW) 

T (M.N) sT(M,N)-(-l.OJO)**NSH*SW ININSH) 

GO  TO  14 

8  I IM,N)=-sWINI2*NSH) 

GO  TO  14 

')  IF  I  (2*NSW)  .EC.  (NW*l) )  CC  TO  1C 
SEP*CA*0FLUAMNW*l-2*NRw) 

T lM,N)x(JCOS<OFlOATlNSHI*THETA))/2,ODO-I-l.ODO)**NSH*SHIN(NSH) 


GO  ro  14 

10  T  (N»N)*(0C0S(CFL0ATI2*NSH1*THETA)  1/2.000 

go  ro  14 

11  SOURC£«O.ODO 
NSELFH=0 

00  12  LDlzl«NH2E 
S  EP*C  A  *QF  LQA  T I  Nk+  1-U)  1-NRMl 
SOURC£*SQURCE+SWIN(NS£LFH  1 

12  CONTINUE 

00  13  L02*l ,NW20 
IF  (L02.E0.NKH)  GO  TO  13 
S£P»CA*CFLOAT I LC2-NRW 1 
SOURCE  *SOURCE*S  WIN! NSELFH) 

13  CONTINUE 

T (M,N) ^SOURCE 

14  CONTINUE 

15  IFILSULVE(NSIZE))  GO  TO  24 

C  CALCULATING  THE  NORMALIZED  RESISTANCE 
RESN= 1.0 
NCS»NW2E*NH 
00  16  NCS  1<  NCS 

RESN«RESN*(  (TINC.NAUo)  )**2J/I:LCAT(NW) 

16  CONTINUE 

IF  (NH20.EG.NW2l)  GO  TO  18 
N0B*NCS*1 

00  17  N0=*N06,NSIZE 

RESN-RESN*! I T (NC.NAUG) >**2I/F LCAT I2*NW) 

17  CONTINUE 

18  WRI Tt ( 6«19)  RESN 

IQ  FORMAT  I /40X23H  NORMALIZED  RESISTANCE  ,F7.4) 

WK1TEI6.20) 

20  FORMAT ( /40X31H-THE  HARMONIC  COEFFICIENTS  ARE-) 

00  23  L*1,NW2C 

M*Nri*( L-l 1*1 
N»NH*L 

WRITM6.21)  L 

21  FORMAT (51X6H  HIRE  .11) 

WRI TE ( 6.22 )  ( T ( I . NAUG) ,  l*M» N) 

22  FORMAT  (5X.12F1C.5) 

23  CONTINUE 
GO  TO  26 

24  WRI  TEI 6 ,25) 

25  FORMAT ( // 53H  THE  T  MATRIX  IS  SINGULAR.  NU  UNIQUE  SOLUTION  EXISTS.) 

26  GO  TC'  2 

27  STOP 
END 

C  THIS  FUNCTION  SUBROUTINE  EVALUATES  THF  ANALYTIC  EXPRESSION  FOR  THE  DEFINITE 
C  INTEGRAL  SWlN. 

DOUBLE  PRECISION  FUNCTION  SHINIIHAR) 

DOUBLE  PRECISION  THET A. SEP, PI . I .C A, PS  I . H. B. C. E ,6 
COMMON  THETA, SEP, PI 

50  IF  (SEP)  51.55,52 

51  PSl»OATANIOSlN(THETA)/ (-2 .000* SEP -DCO SI  THETA) ) I 
GO  TC  53 

52  PSI*PI-OAT AN! OS  INI  THETA)/ (2.000*5  EP+OCOSt THETA))) 

53  SaU SORT (4 .000*1  SEP) **2*l • C00*4 .0C0*SEP*DC0S I  THETA) ) 
A*0.50C*1)CQS(THITA-PSJ*DFIUAT(  IHAR-l) ) 

ll  =  (  1  .000*2.  OOC*  SEP* DCOS I  IliET  A )  )*0CCS  I  PS  I*Df  LO  AT  I IHAR ) ) 

C30. 500*0C0S( THET A*PSI*Cf l.UAT  I  IHAR+l)  I 
IH  (IHAR.EO.O)  CO  lU  54 

SW1  l»l  A*(  S»*2)*t*S*CI/(  I  l.iUUO-  S**2)*l-S)**(  IHAR+1)) 


SO  to  H 

S  SkI»«>SHW.Q}9KI/IS*i  1  ^i»-S«*2>] 
ss  idem 
MB 

c  1ms  react  S3 m  sgmoitimc  sclyes  a  system  or  m  lincm  ejuatiuks  in  n  iinknomns 
c  »  stun:  smssiM  udumtich  kith  caum  riverine. 

Lociui  rwcim  uolveiri 

DAJMHlE  MKISICR  SEr.ri.T.Sufl,TEnr,TOLEft 

onn 

S3*»5.cco 

OO  «1  |>1.1 

CO  *1  >u« 

*1  sai»S8M«aastfii.iii 

lQlE«<IS(ilV3fUIITI1l**2l*l.(a'12 

«n»vi 

MIOrl 

ao  m  mi,mi 
raar*3ACSffiK.Kii 
OO  *2  l<tn.N 

If  iat5tI(I.lll.l£.lEiri  CO  TC  *2 

k2  UMilKE 

If  (IEMP.Lf.IA.HI  CO  TO  T9 
If  I  ITEM .£3. *1  CO  TO  A* 

SO  *3 

tearatiK.it 
rw.iisfiirEKr.it 
*3  riiftxr.iiafEKr 
A*  03  AS  I*KM.ft 

TIS.KI»71 1.XI/TIX.K1 
OO  AS  XM.Ifl 
AS  TC2.JI*f(f  .JT-TI2.KIK71K.J) 

AA  CCKfUU 

ifieusitM.aiiJ.E.ra.H9  co  to  to 

T  IA.A#  *  I~T1  ST.ST  11/7  1  N.Ki 
90  *A 

&**-! 

£10  AT  J-l.l 

LMTl-J 

AT  IU.UIMK.tfU-TK.L!*m.Ull 
as  T».*rii*TW.ftrii/r(K.Ki 
AY  LSiXYEs  JFKLSE. 

£3  TG  71 

70  LS£ti€>.TKlE. 

71  KettfU 
E  W9 


/o! 


C  FORTRAN  IV  PROGRAM  FOR  SOLUTION  EV  TFE  METHOD  OF  LEAST  SQUARES 

C 

C 

C  THli  PROGRAM  USES  THE  METHOD  OF  LEAST  SQUARES  TO  SOLVE  A  SYSTEM  CF  EQUATIONS 
C  FOR  THE  COEFFICIENTS  OF  TR1GCNUMETRIC  SERIES.  THE  SERIES  REPRESENT  THE 
C  NORMALIZED  SURFACE  CURRENT  DENSITY  ON  EACH  MIRE  IN  A  SYSTEM  OF  NH  EQUALLY 
C  SPACED.  PARALLEL.  PERFECTLY  CONOUCT ING  H IKES .  USING  THIS  CURRENT  AN 
C  APPROXIMATE  VALUE  CF  THE  NORMALIZED  HIGH  FREQUENCY  RESISTANCE  OF  THE  SYSTEM 
C  IS  CALCULATED.  THE  NUMBER  OF  HARMONIC  TERNS  USED  TO  OESCRTBE  THE  CURRENT  ON 
C  EACH  MIRE  IS  NH.  THE  RATIO  CA  IS  EQUIVALENT  TO  THE  SPACING  BETHEEN  MIRE 
C  CENTERS  OIVIOEO  BY  THE  MIRE  CIAMETER.  FOR  NM  HIRES  AND  AND  NH  HARMONICS  THE 
C  SUE  OF  THE  MATRIX  TIl.J)  MUST  BE  AT  LEAST  AS  LARGE  AS  T(NUl*hH*l»NWl«NH> 

C  WHERE  NWl-NM/2  FOR  NM  fcVEN  AND  NMl*(NM*l>/2  FDR  NM  ODD. 

LOGICAL  L SOLVE 
EXTERNAL  F1.F2, F3.F4.F5 

COMMON  PI .CA.SEPSIM. SEPMUT, SEPMI M/MATRIX/ T( 2tt  .29) //NM20.NM2E.NRM.J 
i.NSWH. NSMH2.NM 
MR! TE 1 6,1 } 

1  FORMAT  1 1H 1 > 

PI*3. 141593 

2  READ! 5,3)  NM.NH.CA 

3  FORMAT  III, 3X, 12, 3X.F5. 31 
MRI TE I 6,4)NM,CA ,NH 

4  FORMAT I///35XI I.2IH  MIRES,  SPACING  CA»,F6.3,4K,  ,I2,10H  HARMON 

IICS) 

NM20*t 2*NM*l*I-l)**INH*l) 1/4 

NM2E*I2*NM-l*I-lI**NMI/4 

NSIZE*NM2C*NH 

NAUG*NSIZE*l 

00  10  NRM*1,NM20 

SEPSIM*CA*FLOAT INW*l-2*NRM) 

00  10  J*l.NH 
NROM«(NRM-n*NH»J 
00  9  NSW* I.NM20 
SEPMUT*CA*FLOATINSM-NRM) 

SEPMIM*CA*FLUATINW*-1-NSW-NKM) 

00  9  NSMH*1 ,NH 

NSWH2«2*NSMH 

NCOL* I NSM-1 ) *NH*NSMH 

DEL*0.0 

Jl=J 

I FINM20.NE.NM2E .AND.NRH.EQ.NW2C)  Jl*2*Jl 
I  FI J.EC.NSMH)  DEL*Pl/8 .0 
I  FI NRM .EQ.NM20.ANO.NM2E.NE.NM20)  GO  TO  7 
Irt  NKM.EO.NSM)  CO  TO  6 

I  FI NSM.E0.NM20. AN0.NM2E .NE.NM2C)  GO  TO  S 
T INRCM , NCOL )*-ALSQl I J.NSMH, FI ) 

GO  TO  9 

5  T INRCM ,NCOLI=-ALSOl IJ.NSMH2.F2) 

GO  TO  9 

o  T INKOw ,NCOL)  =  CEL-0. 5* ALSO  1 1 J , NSMH  ,F3) 

GO  TO  9 

7  IF! NRM. EQ. NSW)  GO  TO  8 

T INRCM »NCOL ) *-0 .5*ALSQII Jl» NS MF.F4) 

GO  TO  9 

H  T INRCM* NCOL)*DEL 
9  CONTINUE 

T INRCM, NAUG)*ALSQl(Ji.J.F5l 
10  CONTINUE 


/ 


a 


jFu.hot.ve ins uen  on  iu  n 
wkt TE I  6, 1 1) 

U  FORMAT (/40X31H-IHE  HARMONIC  COEFFICIENTS  ARE-) 

DC  14  L  - 1 ,N»20 
M=NH*{ l-1  M-l 
n  =Nh*L 

RRITE(6,12)  L 
12  FORMAT  I 31X6K  WIRE  ,11) 

WKJTE(6,13)  (  T  (  I  «NAUi> ) « l  =  H,  N ) 

1»  FORMAT (5X.12F10.5) 

1 A  CONTINUE 

C  CALCULATING  THE  NGRMALlZfcC  RESISTANCE 
RESN-1 .C 
NCS-=Nw2E*NH 
DO  15  NC= 1, NCS 

RES^KESNM  ( T  INC  *NAUOI  )**21 /UCAT  INW) 

15  CONTINUE 

1F(nw2E.EG.Nw20)  GO  TO  17 
NCS*NCSU 

00  16  N0>NCS,NSI2E 

KESN=RESN*( I T (ND.NAUGI) ♦♦2)/ FLOAT I2*NW) 

16  CONTINUE 

17  wRI TE1 6,18)  RESN 

18  FORMAT  I /A0X2  3H  NORMAL  I  ZEE  RESISTANCE  ,F7.4> 

GO  TC  21 

19  WRI Tt( 6,20) 

20  FORMAT I//53H  THE  T  MATRIX  IS  SINGULAR.  N3  UNIQUE  SOLUTION  EXISTS.) 

21  GO  TO  2 

22  STOP 
END 

C  THIS  FUNCTION  SUBROUTINE  EVALUATES  TEE  DEFINITE  INTEGRALS  WHICH  ARISE  IN  THE 
C  ELEMENTS  CF  THE  MATRIX  T. 

FUNCTION  AL  SCI  I U,NSWH, F ) 

EXTERNAL  FI ,F2, F3.FA.F5 
COMMON  PI 
MP=MAXO(J,NSWk) 

THETAI=0.0 

THE TA2=PI/I2.0*FLOATIMP)) 

T  HE  T  A  3  =  PI-THETA2 
DTHEIA=2.0*THETA2 
ALSUI-C.O 
00  30  IP= 1 , HP 

ALSU  1= ALS0I+GAUSS6I THETA  1, THE TA2.F) 

THE  T A 1  =  THETA2 
30  THE TA2=TFETA2*0 THETA 

ALS0I=ALSCt*GAUSS6ITHETA3,PI,F) 

RETURN 

ENO 

C  FI.  F2,  F  3,  FA  AND  F5  ARE  AUXILIARY  FUNCTION  SUBROUTINES  USED  TO  SIMPLIFY  THE 
C  INTEGRANDS  OF  THE  OEFINITE  INTEGRALS  WHICH  ARE  EVALUATED  NUMERICALLY. 

FUNCTION  FI ( THETA) 

CGHMGN  PI ,CA,SEPSIM, SEPMUT, SfcPPIM  »NW20, NW2E ,NRW, J . NSWH, NSWH2 ,NW 
F 1= ( 0. 5*CCS (FLOAT ( J)* THETA)- (-1.0  )**J*SINTGL( THETA, SEPSIM,J))*I SIN 
liGLI THETA, SEP RUT, NSWH)  ♦ (-1 .0 ) **N SWH*SINTGL ( THETA , SEPMI M, NSWH) ) 

RETURN 

END 

FUNCTION  F2 ( THE T A) 

CuFiMUN  PI  ,CA,  SEP  SIM,  SEPMUT.SE  PRIM,  NW2C.NW2E  ,NRW,  J,  NSWH,  NSWH2  ,NW 
F2=  (0. 5«CCS(  FLOAT!  J  I  *  THET  A)- (  -  1 .0  >**J*S  INTGL1  THE  T  A ,  SEPS  IM ,  J))  *(  SIN 
l  TuL  (  IriE (A,SCPRUT,NSWH2) ) 

kET  urn 


/  d  £ 


END 

FUNCTION  F3 ( THOTA) 

COMMON  PI ,C A , SEPSi M, SEPMOT,SE PPt M ,NH20, NW2E ,NRH, J , NSKH, NS  kH2 ,NH 
F3*(-l.0)**J*SINTGL(TH£TA,S£PSIM, J)*COS (FLOAT! NSHH)*THETA)*(- 1.0)* 
1*NSWH*SINTGL{ THETA, SePSIM,NSWH)*CCS(FLGAT(J)*THETA)-2.0*(-l.0)**(J 
i*NSWH) *SINTGL(THETA, SEPSIM, J)*S!NTCL( THETA, SEPSIM, NSHH) 

RETURN 

END 

FUNCTION  F4(THETA) 

CCMMCN  PI ,CA. SEPSI M,SEPHUT.SEPFiM,NW20,NH2E,NRM,J,NSHH,NSHH2,NW 
F4=C0.S(2.0*FLCAT(JI*THEIA)*(S1NTCL(THETA,  SEPMLT,NSMH)*(-1.0)**NSHH 
1*S1NIGL(THETA,SEPMIM,NSHH)) 

RETURN 

END 

FUNCTION  F51THETA1 

COMMON  PI  ,C A,  SEPSI  H,  SEPMt  T.SEPP1  M.NH2G,  NR 20  ,NRH, J.NSHH, NSMH2 ,NH 

2ERU=0 .0 

NULL*0 

FWERC 

DO  4C  K*1 »NW2E 

S£PMOT*CA*FLOAT  (K-NRWI 

SEPMlM»CA*FLOAT (NH*1-K-NRW) 

1FIK.EC.NRM)  GO  TO  40 

F 5«F5*SINTGL( THETA, SEPMUT, NULL J*S  INTGL! THET A, SEPMIM, NULL) 

40  CONTINUE 

SEPMCT*CA*f LOAT (NH20-NRW) 

I F( NH20.NE.NH2E .ANO.NKM.NE.NM  20)  f 5* I F5*S INTGL ITHETA, SEPMUT .NULL )♦ 

I S IN TGI ( THETA, SEPSI H,NULt ) )*l 0  .5*C0St FLOAT! J )*TH£TA)-(-1.0)**J*SINT 
1 GL (THETA, SEPSIM. J)) 

I F( Nm20.EC.NH2E )  F5*( F5*S INTGLI THETA, SEPSIM, NULL ))*(0.5*C0S (FLOAT ( 

1 J l*THETA)— (- ) .0)** J*SI NTGL1 THETA, SEPSIM, J ) ) 

I F( Nm2C.NE.Nh2E .ANO.NRm. EQ.NH20)  F5=F5* (0.5*CUS(FL0AT( J ) *2.0* THET A 

m 

RETURN 

ENO 

C  THIS  FUNCTION  SUBROUTINE  EVALUATES  THE  ANALYTIC  EXPRESSION  FOR  THE  OEFINITE 
C  INTEGRAL  SINTGL. 

FUNCTION  SINTGL(THETA,SEP,IHARI 
COMMON  PI 

50  IF  (SEP)  51.55,52 

51  P SI=*AT AN(  SI N( THETA) /(-2.0*SEP“COS (THETA))) 

GO  TO  53 

52  PSI  =  PI-AT AN(SlN(THETA)/(2.0 *SEP*COS( THETA ) ) ) 

53  S*SOKT (4. C* (SEP)** 2*1 .0*4 .0*S6P*C0S( THETA) ) 

A=0.5*CCS( THETA-PS I *FLUAT ( 1 HAR-1 ) ) 

B*( 1.0*2. 0* SEP* COS (THETA) )*CUS(PSI*FLOAT( IHAR) ) 

C*0 . 5* CCS (THETA «PSI *FLUAT ( I HAR*1 ) ) 

IF  (IHAR.EO.O)  GO  TO  54 

S 1NT  GL=(A*( S**2 )*B*S*C) / ((l.0-S**2)*t-S)**( IHARMI  > 

GO  TC  55 

54  SINTCL«-(B*S*2.0*C)/(S*( 1.0-S**2I) 

55  RETURN 
ENO 

C  GAUSS*  IS  A  FUNCTION  SUBRUUT INE  WHICH  COMPUTES  AN  APPROXIMATE  VALUE  OF  THE 
C  INTEGRAL  OF  F(X)  OVER  THE  INTERVAL  FROM  X*XL  TO  X*XU.  EVALUATION  IS  DONE  BY 
C  MEANS  OF  A  6-  POINT  GAUSS  QUADRATURE  FORMULA. 

FUNCTION  GAUSS6 (XL, XU, F ) 

A*. 5*( XU*XL  ) 

B=XU-XL 
C*. 466234 8*B 

UAUSS6*. 0656622 5* ( F ( A*C )*F( A— C ) ) 


C=. 3306047*8 

GAUiS6=GAUSS6*. 1803808MF  U*C  ) +F  I  A-C) I 
C  =•  1 193096*8 

u AU bS6=B* I GAUS56*. 2339 570*1 FIA+C)  *FI A-Cl  )  I 

RETURN 

ENO 

C  IH13  FUNCTION  SUEKCUT INC  SOLVcS  A  SYSTEM  OF  N  LINEAK  EQUATIONS  IN  N  UNKNOWNS 
L  BY  US1  Nli  GAUSSIAN  ELIMINATION  WITH  COLUMN  PIVOTING. 

LOGICAL  FUNCTION  LSULVL(N) 

COMMON  /MATRIX/A(2d,29l 
oJ  bUM=C.O 

OL  ot  1=1. N 
OU  o  1  J=  t  »i> 
til  SuM=SUM*ADS< A( 1  ,JI  I 

TUL£k=  <SUN/FLCATINI**2 l*l.JL-6 

NPl  =  N»  1 

NM1 =N-  1 

UO  66  K=l,NMl 

KPl =K ♦  1 

T  EMP  =  A8S( AIK.KJ ) 

I T£MP=K 
00  62  I=KP1 . N 

IF  (ABSIAUtKl)  .LE. TEMPI  GO  TO  62 

T£MP=A8S( A( I ,K)  I 

ITfcMP=I 

62  CONTINUE 

IF  I TEMP.LE.TOLEKl  Gd  TO  70 
IF  ( ITFMP.EQ.KI  GO  TO  6A 
00  63  I=K.NP1 
T  EMP* A  I K »  I) 

A(K.()=A(IT6MP.tl 

63  All TEMP.I >=TEMP 

64  00  65  I=KP1.N 

AH  ,KI=A<  I.KI/AIK.K) 

00  65  J  =  K  PI . NPl 

65  All  «J)=A(  I.J  l-A ( I . K I *A I  X  iJI 

66  CONTINUE 

IFI ABS(ACN.NI). IE. TOLER!  GU  TO  70 

A(N,NP1)=A(N,NP1I/A(N.NI 

00  68  1=1. NM1 

K=N-I 

00  67  J=l.l 
L=NP1- J 

6  7  A(N,NP11  =  A(K,NP1I-A(K,U*A(1,NP1> 

68  AIK,NP11=A(K,NP1I/4IK  K) 

69  L  SOLVE*. FALSE. 

GU  TO  71 

70  L  SOLVE  =  •  T RUE • 

71  RETURN 
ENO 


The  evaluation  of  i'  (o,  tn-1,  p)  closely  follows  that  for  I(o,  m-l,  p) 
carried  out  in  Appendix  A.  Using  standard  trigonometric  identities 
the  numerator  and  denominator  of  the  integrand  are  reqritten  as 


N  =  i  j-A'  cos  [  (p~i)(a'  +  ty)]  +  B'  cos  [  p(o'  +  ty)] 

-  C'  cos  [  ( p+1) (o 1  +  ill)] 

■\ 

-  E'  sin  [  (p-l){9 1  +  iji )]  +  F'  sin  [  p(9'  +  ty)]  - 

G'  sin  [  (p+l)(o'  +  *)]  J 

D  =  s2  +  1  +  2s  cos(9'  +  ty) 

(C-2) 

where 

A'  =  sin(o  -  (p-1)  ijj) 

(C-2a) 

B'  =  -2(1  +  2(m-l)(c./a)cosa)  sin(pi<i  ) 

(C-2b) 

C  =  -sin(9  +  (p+lH) 

(C-2c) 

E'  =  cos  (6  -  (p-1)  i|i) 

(C-2d) 

F'  =  2(1  +  21  cos9)  cos(piji) 

(C-2e) 

G'  =  cos(o  +  (p+1)  ty) 

(C-2f) 

s  =  s2  +  1  +  2s  cos(9'  +  ty) 

(C-3a) 

-106- 


f  n  -  tan 

♦  aV -I, 


-1  /  sini 


(m-i)(cTaV  +  coso  )  iMsl>2"* 
(m-jQ(c/a)  +  cos0  )  ’  =  ~1,  'Z'  "" 


(C-3b) 


The  principle  value  of  tan-i  is  used  in  (C-3b). 

With  equations  (C-2),  (C-3)  and  the  new  variable  (o'  +  \if) 
(C-l)  becomes 


I'  (o,  m-i,,  p)  = 


i=“TT 


-A  cos[  (p-l)$]  -t  B  cos  Ct 
s  +  1  +  2s  cos  (£) 


cos  [  (p-1)  >1  +  E  sin 


sin  [  (p  +  1)$]  do' 


(C-4) 


The  sin  [  (  )  $]  terms  integrate  to  zero  and  the  remaining  terms 

are  in  the  form  cf  a  definite  integral  which  is  evaluated  in  Appendix  . 

A,  equation  (A -7).  With  these  integrations  performed  I'  (q,  m-i,  p)  is 


I'  (a,  m -i,  p)  = 


(1-  «  K-s)1 


[  A's  +  B's  +  C'  ]  ,  p  =  l,  2  -- 


(C-5) 


Appendix  D 

Listing  of  Computer  Program 

This  appendix  contains  a  computer  program  written  in 

Fortran  IV  language  for  use  on  the  I.  B.  M.  360/65  computer. 

The  program  computes  the  coefficients  a  and  b  of 

mcp  mcp 

the  trigonometric  series  for  the  normalized  correction  current 
densities  g  . 


/ 


f 
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C  FORTRAN  IV  PROGRAM  FOR  CALCULATING  THE  CURRENT  DISTRIBUTIONS  CAUSEO 
C  BY  INTERCONNECTING  HIRES  IN  THE  EXPERIMENT AL  MODEL* 

C 

C 

C  THIS  PROGRAM  USES  THE  METHOD  OF  COLLOCATION  T3  SOLVE  A  SYSTEM  OF 
C  EQUATIONS  FDR  THE  COEFFICIENTS  OF  TRIGONOMETRIC  SERIES.  THE  SERIES 
C  REPRESENT  THE  NORMALIZED  SURFACE  CURRENT  DENSITY  ON  EACH  WIRE  IN 
C  SYSTEM  OF  NW  HIRES.  THE  NUMBER  OF  HARMONIC  TERMS  USED  TO  DESCRIBE 
C  THE  CURRENT  UN  EACH  WIRE  IS  NM.  THE  RATIO  CA  IS  EQUIVALENT  TO  THE 
C  SPACING  BETWEEN  WIRE  CENTERS  DIVIDED  BY  THE  WIRE  DIAMETER.  THE 
C  NORMALIZED  DIMENSIONS  C AH 2  AND  CHL2  ARE  THE  TUBE  RADIUS  DIVIDED  BY 
C  THE  TUBE  HALF  HEIGHT  AND  THE  TUBE  HALF  HEIGHT  OIVIDED  BY  THE  LENGTH 
C  OF  THE  INTERCONNECTING  WIRES. 

LOGICAL  LINECN 

COMMON  THETA.SEP.PI/MArRIX/A(36t37) 

WRI  TEI6.1 ) 

1  FORMAT  (HI) 

PI*3. 141593 
CAH2«0.048645 
CHL  2=0.3.0824 

2  READIS.3)  NW.NH.CA 

3  FORMAT (U.3X.12.3X.F5.3) 

WRITE(6,4)  NW.CA.NH 

4  FORMAT I///35XI1 ,21H  WIRES.  SPACING  CA-.FS.3.4H.  .12.10H  HARMON 

1ICS) 

NW21«(2*NH+U(-l)*«(NW*l))/4 
NW2  =  l2*NW-lH~l>**NW)/4 
NSI ZE=NW21*NH 
NAUG=NH*iNH21+l 
00  17  L-l.NH 
NR1= (NW21-1 ) *NH+L 
C  SETTING  COLLOCATION  POINTS 

THETA=PI*2.0*<FL0AT(L))/IFI  QATCNH«-1)> 

DO  17  M=U.NR1.NH 

NRH=1.IM*L)/NH 

IF  (NW21.EQ.NW2)  GO  TO  5 

I  FI NRW.ED.NW21 >  THETA-I THETA-PI ) /2.0 

IF( THETA. LE. 0.0)  THET A=2.0*P I .THETA 

5  NSW=1 
N1*NH 

DO  17  N=l  ,NAUG 
IF  (N.EO.NAUG)  GO  TO  15 
IF  (N-Nl)  7.7,6 

6  Nl-NUNH 
NSH-NSW.l 

7  NSH-N-NHMNSW-1 1 
IF(N.EQ.2*(N/2) )  NSH«NSH/2 

IF!  N.NE.2  *(N/2) )  NSH« <NSH*i  )/2 
IF  ( NRW.cQ. NSW)  GO  TO  11 
SEP=CA*FlOAT«NSH-NRH) 

IF  ( (2*NS W) • EQ. (NW+l ) )  GO  TO  9 
I  F(  2*( N/2 ) . EQ.N)  GO  TO  8 
AIM.Nl— AINEINSH) 

SEP-CAOFLUAT (NW.l-NSW-NRW) 

A(N »N>  =A(  M»N)-(  ”1.0)**NSH*AINE (NSH) 

GO  TO  17 

8  A(M.N) ■-AINES(NSH) 

SEP«CA#FL0AT(NW.1-NSW-NRN) 

A(M,N)-A(M,N).(-1.0)**NSH*A1NES(NSH) 


■i 

a 

i 

-j 

i 
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GO  TO  17 

9  I  FI  2*1  N/2).  EQ.  lt  GO  TO  10 
AIN.N) — AINE 12+NSH) 

GO  TO  17 

10  AIM.N)— A1NESI2*NSH-1> 

GO  TO  17 

11  IF  II2*NSM).E0.  INW+ll)  GO  TO  13 
SEP«CA*FL0ATINW+l-2-NRW) 

IF(2*IN/2).EQ.N>  GO  TO  12 

AIN. N)*ICOSI  FLOAT!  NSH)*THET4)  1/2.0-1-1.01  **NSH*AINE INSH) 

GO  TO  17 

12  AIM .N! -( SIN! FLOAT) NSH) -THETA)  ) /2. 0+1-1. 0) **NSH*AINESINSH) 

GO  TO  17 

13  1 F!  2*1  N/2 1.E0.N1  GO  TO  14 

AIN, N) -1C JSI FLOAT! 2-YSH) -THETA) 1/2.0 
GO  TO  17 

14  A(N«N)-($iNI  FLOAT!  2 -N Shi- H* THETA)  1/2.0 
(»0  TO  17 

15  SOURCE -0.0 

00  16  LOl-l.NH 
$EP-CA*FLdATIL01-NRW) 

S0URC1  -  ICAH2/SQRT! 1.3+ ICHL2 )**2+4<  0*! SEP*CAH2*CHL2 )**2t)*t!C 
iHL2 )**2*l SIN! THETAI “2.0*SEP*CAH2*CHL2*C  OS! THETA) )/ !  i.0+4.0-1  SEP*CA 
i H2*CHL2)**2 )+Sl NITHET A)  /  (1»0+4«0*ISEP*C  Art2)  **2 )  1-1 1 l.0+2.0*$EP*C0S 
II  THETA)  1/1 1.0+4.0*S£P**2+4, 0*SEP*C0S!THETA) 1 1*11. 0-1. O/SQRT!  1.0+4. 
10*ICAH2)**2*SEP**2) 1 
SOURCE' SO JRCE+SOURCl 

16  CONTINUE 
AIM, NI-SJURCE 

17  CONTINUE 

18  IF  IL1NE3NINSI2EII  GO  TO  23 
WRITEI6.19) 

19  FORMAT  I /+0X31H" THE  HARMONIC  COEFFICIENTS  ARE-) 

DO  22  L-1.NWZ1 

N«NH*I L-l t  +  i 
N-NH-L 

MRI TEI 6. 20)  L 

20  FORMAT { 51X6H  WIRE  ,11) 

WRITE!6,21I!AII,NAUG),  l-H.N) 

21  FORMAT  1 5X, 12 Flo. 5) 

22  CONTINUE 
GO  TO  25 

23  WRI TEI 6,2+) 

24  FORMAT !// 53H  THE  A  MATRIX  IS  SINGULAR.  N3  UNIQUE  SOLUTION  EXISTS.) 

25  GO  TO  2 

26  STOP 
END 

C  THIS  FUNCTUN  SUBROUTINE  EVALJATES  THE  ANALYTIC  EXPRESSION  FOR  THE 
C  DEFINITE  INTEGRAL  A1NE. 

FUNCTION  AINEIIHAR) 

COMMON  THETA, SEP, PI /MATRIX/ A I  36,37) 

50  I  FI  SEP )  51,54,52 

51  I  FI  THETA.  LE.  PI  I  PSI  — ATAWISIN ITHETA)/ 1-2, O-SEP-COS  ITHETA) ) ) 

IF! THETA. G£. PI)  PSl- ATANI -SIN ITHETA) /!* 2, 3+SEP-COS! THETA) ) ) 

GO  TO  53 

52  1FI  THETA.LE.PI)  PS  I-P I+ATANI  S  INI  THETA)/ 1 2.0-SEP+COSITHETA) )) 

IF! THETA. GE.PI)  PSl -PI-ATAN1 -S  IN  ( THETA) /!  2, 0+ SEP+COS (THETA) ) 1 

53  H-ll. 0+2, 0*SEP*COS!  THETA)  )*CUS!PSI*FLJ»Ti  IHAR  t ) 

B— 0.5*CJSI THETA-PSI*FLOATI IHAR+1 ) ) 

C—0.5*CJS<  THETA+PS1  -FLOAT!  IHAR-l ) ) 

E»4.0*l  SEP)  **2+ 2*0+4.  0*SEP*CUSI  THETA) 


F«I2.0*Ie-1.0)«*0.5)/£ 

G»H1.0-F**2)**0.5-1.0)/F 

AINE«tG**tiHAR-ll)*<H*GHB*tG**2)*C)/IE*U1.0-F**2>**0.5>) 

54  RETURN 
ENO 

C  THIS  FUNCTION  SUBROUTINE  EVALUATES  TFE  ANALYTIC  EXPRESSION  FOR  THE 
C  DEFINITE  INTEGRAL  AINES. 

FUNCTION  AINESI 1HAR) 

COMMON  TH£TA,$EP,PI/NATR1X/AI36,37) 

80  I  FI  SEP  I  8 1 , 84 , 82 

81  IFITHETA.LE.PII  PSI=-AT ANISIN!  THETA)/ 1-2.  Q*jEP-COSt THETA) ) ) 

IF!  THETA. GE. PI)  PS MATANi-SINI THETAI/I-2. (WSEP-COSt THETA)  >  ) 

GO  TO  83 

82  IFITHETA.LE.PII  PSI=PI*ATANI SINI THETA)/ 1  2. 3*S EP*COS ITHETA)  I) 

I  FI  THETA. OS. P 1 1  PS  I =P I-ATANI - S INI  THETA) /1 2. 0*  SEP+CQS I  THETA ) ) ) 

83  H1=I1.0»2.0*S£P*COSITHETA))*SINIPSI*FLOATUHAR>) 

8 l=-0. 5*5 INI PSI *FL3AT 1 1 HAR+1 l-THETA) 

Cl=-0. 5*5  INI PSI*FL3AT 1 1HAR-1 ) HTHETA) 

£=4.  0*  l  Si  P)  **2+2. 0*4.  0*SEP*C0St  THETA) 

F=I2.0*(c”1.0i**0.5)/E 
G*I  1 1. 0-F**2 1**0.  5-1. 0)/F 

AIN£S*IG**I 1HAR-1) )*tHl*G*Sl *1 G**2  MCI  1/ 1 E*I I l.Q-F**2)**0. 5) ) 

84  RETURN 
cNO 

C  THIS  FUNCTION  SUBROUTINE  SOLVES  A  SYSTEM  OF  N  LINEAR  EQUATIONS  IN 
C  N  JNKNOWNS  dY  USING  GAUSSIAN  ELIMINATION  WITH  COLUMN  PIVOTING. 
LOGICAL  FUNCTION  L1NEQNIN) 

CCKMON  THETA.SEP.PI/KATR1X/A I  36, 37) 
oO  SUM*0.  0 

DO  61  1=1 .N 
00  61  J=l.N 

61  SUN*SUMHAdS<AII,J)l 

i UkER> I $UM/FL0ATlNI**2)*l.0E-6 

NPl«N+l 

NM1=N-1 

00  66  K=1 .NM1 

KPl=K+l 

TEMP»ABSI  AIK.lO  I 

ITEMP-K 

00  62  I-KPl.N 

IF  lABSUIl.KD.LE.TcMP)  GO  TQ  62 

T  EMP»ABSI  AII.K)  I 

ITEMP=1 

62  CONTINUE 

IF  I  TEMP.  LE. TOLER)  GO  TO  70 
IF  I ITEKP.EQ.KI  GO  TO  64 
00  63  I »< , NP1 
TEMP=A (K. 1 ) 

A  IK ,  I )  =  AI  1TEMP.  1 1 
83  A 1 1  TEMP, l  MTEMP 

64  00  65  !«*P1.N 

All  «K)  *AI  I.KI/AIK.K) 

00  65  J=RP1 , NP1 

65  All ,JI*AII,J)-Atl,K)*AIK,J) 

66  CONTINUE 

IFIABSIAIN.NM.LE. TOLER)  GJ  TO  70 
AIN.NPIMAIN.NPU/AIN.N) 

00  68  1=1 »NM1 
K-N-i 

00  67  J«i,I 
L*NP1-J 


67  AIK.NP1MA(K,NP1)-A<X,L)*AIL,NPI> 

68  AIK»NP1I*A(K,NP1)/AIK»KI 

69  L  INECN*. FALSE. 

GO  TO  71 

70  L INEUN-.TRUE. 

71  RETURN 
ENO 


/  / ' 
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